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SOMMAIRE 

On a tenu un atelier les 9 et 10 juillet 1987 dont l'objectif 
etait d'elaborer un sous-modele de la chaine alimentaire 
aquatique qui serait incorpore au modele du sort des polluants et 
au modele du transport des polluants qui font partie de la SMID. 
Au cours des discussions en groupe et des exposes faits par les 
conferenciers invites, on a degage un certain nombre de concepts. 
En ce qui concerne la model isat ion de la chaine alimentaire, les 
produits chimiques qui s •accumuleraient dans les organismes 
vivants ont ete isoles en tant qu' elements critiques. 

Les participants aux seances en petits groupes ont convenu de 
creer un modele de base de la chaine alimentaire, qui 
comprendrait quatre ou cinq niveaux trophiques bases sur la 
taille des organismes. lis ont egalement convenu d'utiliser ce 
modele comme instrument de tri des produits chimiques pour 
determiner s'il faut se preoccuper, dans des endroits precis, du 
processus de bio-accumulation. Les participants ont insiste sur 
1* importance de commencer 1 'elaboration d'un modele specif ique 
pour un site controle, a partir des donnees connues. 



OPENING REMARKS 



YOUSRY HAMDY 
ONTARIO MINISTRY OF THE ENVIRONMENT 
WATER RESOURCES BRANCH 
TORONTO, ONTARIO, CANADA 
M4V 1P5 



Welcome to the Food Chain Modelling Workshop. 

In June 1986 the Ontario Ministry of the Environment announced the 
Municipal Industrial Strategy for Abatement (MISA) program in order to 
eliminate discharges of persistent toxic contaminants to our waterways. 
The concept of virtual elimination stated in the MISA white paper focused 
on the challenge of dealing effectively with the increasing presence of 
toxic contaminants in Ontario's lakes and rivers. 

Events such as those associated with the perchloroethylene spill in 
the St. Clair River and the presence of dioxin and mirex in the Niagara 
River prompted the attention of scientists to the fate and behaviour of 
toxic contaminants in the aquatic environment. 

In this workshop we will focus on the fate of organics such as 
chlorinated benzenes, PCBs, polyaromatic hydrocarbons and dioxins. In 
general, such organic toxic contaminants pose a risk to fish, plants and 
other wild life and even man at low concentrations. They are persistent, 
accumulate in the environment and can be biomagnified in the food chain, 
thereby resulting in fish contaminant levels which may be thousands of 
times higher than in the water. 

The development and application of food chain modelling serves as a 
valuable tool to provide a capability for testing hypotheses, identifying 
missing or poorly understood processes and integrating factors related to 
the fate of toxics. These models will provide service throughout the 
entire resource management process beginning with the problem definition 
and culminating in the selection of effective management options. 

As a minimum requirement, the MISA white paper indicated that the Best 
Available Technology Economically Achievable (BAT-EA) will be used to 
establish effluent limits for industrial and municipal sectors. The 



application of BAT will not always provide sufficient control of 
discharges to prevent adverse effects on the ecosystem. 

Ve therefore must use models as tools to predict with a high degree of 
accuracy the fate and effect of loadings and concentrations of trace 
organics on each ecosystem compartment and to ensure compliance with water 
quality objectives and guidelines for aquatic life, wild life and public 
health over the short and long term. Monitoring of ecosystem compartments 
will require the collection of data on effluent quality, water quality, 
sediments (both active and non-active layers), phyto- and zoo-plankton, 
benthic community (taxa and body burden) and fish. Ve hope that the total 
picture can be obtained from an in-depth understanding and description of 
the food chain and the associated model application. These descriptions 
are the goal of this workshop. 
ORGANIZATION OF THE WORKSHOP 

In order to focus the workshop, several formal seminars were presented 
on the first day of the workshop (see attached programme). These seminars 
addressed both the fundamental questions of uptake as well as 
state-of-the-art food web modelling. The formal seminars were followed by 
work sessions. These were organized into three groups, namely: 

1) plankton 

2) fisheries 

3) sediment/benthos. 

The participants and coordinators of these work groups are shown on 
the attached list. Dr. McCorquodale and I participated in all three 
groups. 

A coordinators meeting was held at the end of the first day. 






The coordinators reported the findings of their respective group to 
the groups as a whole on the final day of the workshop. These findings 
and the ensuing discussion are summarized by J. A. McCorquodale in Chapter 
8. 

Several handouts (see samples attached in Appendix A) were prepared to 
help focus the discussion of the work groups. 
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MISA WATER QUALITY TRACK 



C.F. SCHENK 
WATER RESOURCES BRANCH 
ONTARIO MINISTRY OF THE ENVIRONMENT 
TORONTO, ONTARIO, CANADA 
M4V 1P5 



INTRODUCTION 

In June 1986 the Ontario Ministry of the Environment announced the 

MISA - "Municipal-Industrial Strategy for Abatement" - program. This new 
approach to controlling point-source water pollution has an ultimate goal 
the virtual elimination of toxic contaminants from all industrial and 
municipal effluents discharged into the province's waterways. The program 
is being developed in consultation with Environment Canada, industries, 
interest groups and the general public. 

For the first time in Ontario, there will be regulations to limit 
loadings of trace organic and inorganic substances to surface waters. For 
each of the major sectors of industry and the municipal sector, two 
regulations will be developed. Initially, dischargers will be subject to 
Monitoring and Reporting Regulations which will require them to identify, 
measure and, report on concentrations and amounts of toxics in their 
effluents. Minimum legal requirements for the submission, accuracy and 
reliability of self-monitoring information, including specifications for 
sampling and analytical protocols and quality assurance and quality 
control procedures, will be set. This self-monitoring program will be 
audited by the Ministry of the Environment. 

This initial monitoring phase will be followed by Compliance or Limits 
Regulations that will specify allowable concentrations, as well as amounts 
or loadings, of toxic pollutants for each discharger. This will be done 
on a sector-by-sector basis and the specified levels will be based on the 
Best Available Technology Economically Achievable (BATEA) . 

Joint technical committees for each industrial sector and the 
municipal sector will negotiate practical and effective requirements for 
each regulation. 



More stringent effluent limits may be required for specific sites to 
protect sensitive water bodies. This action would occur in cases where 
water quality assessments show that controls beyond those required by 
BATEA may be necessary to protect the environment (Fig. 1.1). 

The BATEA and Water Quality tracks are complementary in nature. This 
interrelationship will enable the implementation of effective controls 
with impact on the environment being the final determinant. 
MISA WATER QUALITY TRACK 

The Ministry is undertaking a major review and revision of its 
traditional water quality impact approach to enhance the province's 
capability to assess ecosystem response to municipal and industrial 
discharges. 

These evaluations include: 

i) The development of ambient water quality objectives (PWQO) for 
many additional toxic substances. This will involve updating 
"Blue Book" ambient numerical and narrative requirements to 
protect water quality; (Fig. 1.2); 

ii) The revision of policies and procedures for the Water Quality 
based approach to pollution control; 

iii) The development of testing techniques to assess the impacts of 
toxics at selected pilot sites; 

iv) The development of precepts for controlling toxic discharges; 

v) The definition and staging of virtual elimination of toxic 
contaminants entering our waterways; 

vi) The development of a relatively simple screening process to 
determine whether BATEA-controlled discharges within a given 
sector will meet receiving water requirements; 

vii) The selection and validation of biological effluent testing 
procedures to determine acute and chronic toxicity, genotoxicity 
and bioconcentration; 

viii) The development of acceptable guidelines for sediment and biota. 



The competition of these components will allow identification of 
situations where more stringent effluent limits are called for, the 
application of approval assessment procedures and development of effluent 
requirements based on water quality impacts on a site-specific basis. 

Following an initial screening process, discharges requiring more 
in-depth receiving water assessments will be subject to intensive 
site-specific modelling efforts to establish local impacts and related 
controls. Fate and transport models now evolving through exercises such 
as this workshop will be essential to determine far- field concentrations, 
partitioning and potential aquatic effects. 

Looking to the future, whole lake effect modelling must be integrated 
with toxocity and bioconcentration assessments as a basis for 
prioritization of sources and related controls. Satisfactory and points 
for toxic contaminants must be established to permit trend- through- time 
evaluations at critical locations as a means of monitoring the success of 
remedial efforts and maintaining vigilance on the integrity of the aquatic 
ecosystem. 
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FIGURE CAPTIONS 



Fig. 1.1. Technology and water quality tracks. 



Fig. 1.2. The revised blue book will be the expression of the Ministry's 
Water Quality Management Program including MISA's principles, 
policies, objectives, and ambient and discharge requirements. 
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FIG. 14 Technology and Water Quality Tracks 



The Water Quality Track is one of 
the two major components of MISA 
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FIG. 1-2 The revised Blue Book will be the expression of the Ministry's 
Water Quality Management Program including MISA's principles, 
policies, objectives, and ambient and discharge requirements. 



CHAPTER 2 

A MODELLING FRAMEWORK FOR ESTIMATING 
CONTAMINANTS LEVELS IN BIOTA 



J. A. McCORQUODALE 

DEPARTMENT OF CIVIL ENGINEERING 

UNIVERSITY OF WINDSOR 

WINDSOR, ONTARIO, CANADA 

N9B 3P4 



OBJECTIVE 

This workshop was organized to exchange information between Great 
Lakes aquatic biologists and water quality modellers. The objective of 
the workshop was to suggest a food chain/web modelling framework that 
could be developed and tested within the MISA programme. 

Figure 2.1 shows a sample of a generic food web of the type that can 
be handled by state-of-the-art "food chain" models. Figure 2.2 shows a 
chart for the preparation of a predator-prey diet distribution matrix. 
Specific questions that were directed to each work group are shown in 
Appendix A) . 
GENERAL BACKGROUND 

The early water quality models were developed to answer questions 
related to the definition of zones of severe or even toxic impacts. These 
models considered primarily the near field mixing and far field dispersion 
processes. Questions of partitioning to solids and uptake in biota were 
not considered. Such models were and still are useful in evaluating 
Limited Use Zones (LUZ) for conventional pollutants, e.g. TDS, un-ionized 
ammonia, and residual chlorine. However, these models are completely 
inadequate for organic toxics and certain heavy metals which partition to 
the sediment and are accumulated in the biota. 

A new generation of models are now being developed which eventually 
will allow us to trace the fate of a contaminant from its source to the 
top of the food chain. These new models require a vast data base which 
includes measurements of mass loading and concentrations in ambient water, 
suspended sediment, bed sediment as well as in the biota at all trophic 
levels. In addition, the physical and chemical properties of the 
contaminant must be known. 



MODELLING FRAMEWORK 

The first step in any modelling exercise is the selection of the 
appropriate model that will best answer the questions related to a 
particular management decision; for example, what is the load allocation 
which will satisfy the provincial water quality objectives for ambient 
water, sediment, and sports fish? 

All higher level surface water quality models consist of three 
components: 

a) an hydrodynamic submodel, 

b) an exposure submodel, 

c) a biota uptake submodel. 

A systems diagram showing these components is given in Fig. 2. 3. 
HYDRODYNAMIC SUBMODEL 

Only a brief description of hydrodynamic modelling can be included 
here. These models may be one, two- or three-dimensional in space and may 
be steady or unsteady. There are a number of available and documented 
models that can be used in the Great Lakes System. The US EPA supports a 
general unsteady network type hydrodynamic model called DYNAHYD3 (Ambrose 
et al. 1986). A simple two-dimensional steady-state model for large 
rivers has been developed at the University of Windsor (McCorquodale et 
al. 1987); this model is referred to as the k-E model. It is based on the 
work of Rastogil and Rodi (1978) and Krishnappan and Lau (1983) and 
provides velocities and dispersion coefficients on a variable grid across 
and along a river. There are a number of finite element and finite 
difference models available for computing circulation in the Great Lakes 
(Schwab et al. 1981, Simmons 1971 and Ibrahim and McCorquodale 1985). 
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Exposure Submodels 

Exposure models are mass balance models which attempt to account for 
transport, partitioning, degradation and loss of any contaminant 
introduced to the system (see Fig. 2.4). These models are meant to give 
estimates of the ambient contaminant levels in all compartments of the 
lake or river. 

In order to accomplish this, it is essential to know the physical and 
chemical properties of the pollutant. Some of these properties are: the 
octanol-water partition coefficient, the organic carbon partition 
coefficient, the molecular weight of the chemical, the Henry's law 
constant, the vapour pressure of the compound, the specific gravity, and 
the aqueous solubility of the toxicant chemical species. Other useful 
information for prediction are: bacterial population density in water, 
bacterial population density in bed sediments, concentration of solids in 
bed sediments, concentration of solids in the water column, the size 
distribution of the solids in the water column and the size distribution 
of the solids in the river bed, and the organic carbon content of the bed 
sediment. 

Generally speaking the volatilization and decay mechanisms are 
important for long term considerations especially when flow enters slow 
moving bodies such as lakes. A thorough discussion of decay mechanisms is 
given in the EPA Reports dated June 1985 and March 1983. The major 
mechanisms include: volatilization, photolysis, hydrolysis, 
biodegradation, oxidation, and reduction. 

There are several documented exposure models that have been applied 
for the Great Lakes System. Some of these are shown in Table 1. A 
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typical example is the EPA TOXIVASP model (Ambrose et al. 1983) which 
incorporates the features of WASP (DiToro et al. 1981) and EXAMS (Burns et 
al. 1982) including sediment-water interaction. Moreover it had already 
been used for toxic chemicals in a Great Lakes environment (Great Lakes 
Institute 1984). 

The hydrodynamics for TOXIVASP can be obtained from DYNAHYD3 or from 
the k-E model described above. TOXIVASP includes the main mass transfer 
mechanisms to account for the interactions between dissolved chemical both 
with suspended solids and with stationary sediments in the river bed. 
These mechanisms, as shown in Fig. 2.5, include: 

a) advection and dispersion of dissolved contaminants, 

b) sediment contributions from point/nonpoint sources of solids and 
contaminants into the water, 

c) volatilization and biological degradation. 

The TOXIVASP model treats volatilization by the two film concept that 
is used by its predecessor, EXAMS (Burns et al. 1982). 
The TOXIVASP model is a cell model which assumes that: 

a) all segments are well mixed, 

b) sorption is an instantaneous process within each segment, 

c) the chemical properties of the compound can be coupled with the 
characteristics of the environment to formulate a first order 
reaction for the degradation process, 

d) all the first order rates can be combined linearly. 

Vithin each cell the elementary mass balance equation is solved to 
yield the time variable mass and concentration in each cell. Figure 2.6 
shows schematically the partitioning of the total contaminant between the 
dissolved and particulate phases. 
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Figure 2.7 shows a typical schematic of the TOXIVASP cell arrangement: 
the number refers to the water column cells, the corresponding active and 
deep bed cells. 
The Food Chain Submodel 

The purpose of the food chain model is to relate the body burden in 
various trophic levels within the food chain to the prevailing or 
anticipated ambient contaminant levels. 

Food chain models vary greatly in sophistication and data 
requirements. The following are some of the choices to be made in the 
selection and development of a model to predict body burdens in biota: 

a) simple partitioning (BCF x ambient concentration), or a food web 
including uptake from water and food sources, 

b) the model can be either steady or unsteady (dynamic), 

c) in addition, a particular species can be treated by age or size 
classification, 

d) the food chain submodel can be uncoupled or coupled to the 
exposure model, 

e) the food chain model can be treated as a pathway (predator/prey) 
system or as a complete energy and mass balance system. 

One of the simplest biotic contaminant models was proposed by Whitmore 
(1977) as illustrated in Fig. 2.8. This model uses simple partitioning 
between the dissolved and biotic phases. 

Neely (1977) proposed a similar model (illustrated in Fig. 2.9) which 
attempted to predict contaminant levels in fish based on uptake and 
clearance rates. His model included a reasonable representation of the 
fate of a contaminant in the water phase. 

Chapra (1982) proposed a comprehensive food web model as illustrated 
in Fig. 2.10. This model is generic in nature but includes both pelagic 
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and benthic components of the food chain. He included the following 
transfer processes: 

a.) flushing, 

b) loading, 

cj settling, 

dj diffusion, 

e) burial, 

f\ harvest, 

g) sorption, 

n) feeding, 

i) uptake, 

j) release, 

k) egestion, 

1) vaporization. 

Manhattan College (Connolly et al. 1986) has developed a comprehensive 
modelling framework for a food chain which is represented schematically in 
Fig. 2.11. Their model (WASTOX) can be run in several modes, for example, 
as a food chain model with steady or unsteady exposure levels obtained 
from field measurements or their exposure model or any other appropriate 
exposure model such as TOXIWASP. VASTOX can accommodate several trophic 
levels (predator/prey relationships) with interaction at a number of 
levels, age classification within a species, migration of species within a 
cell structure, growth of the species, as well as unsteady loads and 
exposure levels. Furthermore, this model can be developed as a generic 
model or with sufficient site data it can be developed as a site-species 
specific model. 

Thomann (1987) gives the rate equation used in the WASTOX food chain 
model as, 



di/. 

- = K -C n + S Q;:C-.U. - K.U; (1) 
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in which v. = biota contaminant concentration in species i 

c D = dissolved concentration 

K ui = uptake or BCF 

a = assimilation efficiency of contaminant, 

C- ■ = consumption of species j by species i 

K. = net growth + excretion rates for species i 
A PROPOSED PATHWAY FOOD MODEL 

In order to progress towards a realistic food chain model the food web 
system shown in Fig. 2.1 is proposed as a starting point for discussion. 
This system can be accommodated by WASTOX provided that we have sufficient 
knowledge about the following: 

a) the exposure levels for the selected contaminant, 

b) the KOV of the contaminant, 

c) the growth characteristics of the species of interest, 

d) the respiration rates for the selected species. 

e) the excretion rates for the selected species, 

f) the lipid fractions of the selected species, 

g) the migration patterns for all selected species, 

h) the predator/prey relationships for all selected species, 
i) biota mass by species and age classification. 

The feeding relationships can be summarized on a consumer food matrix 
chart such as the one illustrated in Fig. 2.2. 

Figure 2.12 lists a number of factors that a food web model should 
include. Many of these need further quantification through field and 
laboratories research before they can be adequately included in a 
mathematical model. 
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FIGURE CAPTIONS 

Fig. 2.1. System diagram for a generic aquatic food chain model. 

Fig. 2.2. Predator-prey diet matrix. 

Fig. 2.3. A systems diagram for a complete surface water quality model 
(Hague, 1980). 

Fig. 2.4. Typical transfer processes in exposure models. 

Fig. 2.5. Processes in the TOXIVASP exposure model (Ambrose, 1983). 

Fig. 2.6. The partitioning process in an exposure model. 

Fig. 2.7. A typical planwise cell discretization for TOXIWASP. 

Fig. 2.8. The Whitmore model (Chapra, 1983). 

Fig. 2.9. The Neeley model (Chapra, 1983). 

Fig. 2.10. The Chapra model (1983). 

Fig. 2.11. The Manhattan College model, WASTOX (Connelly, 1986). 

Fig. 2.12. Some factors that affect the contaminant uptake and 
accumulation in biota. 
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FIG. 2-1 System Diagram for a Generic Aquatic Food Chain Model 
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FIG- 2-2 Predator-Prey Diet Matrix 
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FIG. 2-3 A Systems Diagram for a Complete Surface Water Model 
(Hague, 1980) 
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FIG- 2-6 The. Partitioning Process in an Exposure Model 



FIG. 2.7 A Typical Planwise Cell Discretization for TOXIWASP 
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FIG. 2-9 The Needly Model (Chapra, 1983) 
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INTRODUCTION 
Background 

The work described in this summary is part of an overall research 
project on the variability of chemicals in aquatic food chains (Thomann 
1987a). The specific details of the relationship of food chain behavior 
and the log K of organic chemicals is given in (Thomann 1987b). 

The purpose of this summary is to provide a brief review of the basic 
structure of the food chain model, indicate the relationship to the 
octanol water partition coefficient and draw some general conclusions 
regarding the need for and use of food chain modeling in water quality 
monitoring and allocation activities. 

There has been a considerable amount of discussion in recent years 
about the actual occurrence of food chain accumulation. That is, a simple 
lipid partitioning of the chemical into the lipid pools of aquatic 
organisms has been used to determine the concentration of the chemical. 
This simple bio-concentration calculation has been widely used in the 
United States in establishing water quality criteria for organic 

chemicals. 

In this view, the aquatic system is seen as a series of "lipid 
packets" that sorb the organic chemical directly from the water. Food 
chain effects are non-existent in this model since there is no apparent 
mechanism for chemicals to accumulate from one level of the food chain to 
the next. The procedure is then very simple. All one needs to know is 
the octanol water partition coefficient (K QW ) and the available chemical 
concentration in the water. The concentration in the aquatic organism is 
then simply given by 
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t - c K ou p (i) 

where 

v - organism chemical concentration (/ig/g(w)) 

c = chemical concentration in the water (/ig/L) 

K a the octanol water partition coefficient, and 
ow 

p = lipid fraction of organism (kg(lip)/kg(w)) . 

The difficulty with this simple procedure is that there is some 
considerable growing evidence that for certain organic chemicals, there is 
an apparent significant increase of the chemical concentration in the 
organism above that which would be expected from uptake from the water 
only. 

Since there are data that indicate some significant bioaccumulation in 
the food chain, it is necessary to construct modeling frameworks that 
incorporate the principal features of chemical uptake and food chain 
transfer. A model framework for this purpose is discussed in the next 
section. 

TOE FOOD CHAIN MODEL 
Generic Versus Site Specific Food Chains 

There are two approaches to structuring a food chain model, each with 
differing purposes and objectives. The first is a generic approach where 
generalized relationships are applied to a simple "one-dimensional" food 
chain. That is, a single predator feeds on a single prey. This generic 
type structure is useful for general screening analyses of chemicals and 
general behaviour of chemicals in abstract but very easily understood food 
chains. This is the modeling focus presented in this paper. 

The second approach is to construct a food chain model that is fully 
site specific. The food web is described as completely as 
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possible, using for example, analyses of food contents to infer feeding 
patterns. Organism specific growth rates are used, associated food 
assimilation efficiencies are assigned and any specific migratory behavior 
is explicitly included (Thomann et al. 1984, Connolly et al. 1985). This 
approach is most useful for describing the particular behavior of an 
individual food chain for use in detailed water quality analyses and 
determination of site specific waste load allocations. 
Equations for Generic Food Chain-Water Uptake Only 

Considering a simple food chain of four levels with level #1 as the 
base of the food chain represented by the phytoplankton. Then a mass 
balance of a chemical around a defined food chain compartment is given by 
(Thomann 1987) 

-rr 1 = k'. w. c + a. ■ 1 C^ • 1 w- -±± - iy, (2) 

dt ui l 1,1-1 i,i-l i w j_l 



i i 



where 

v'- - chemical whole body burden (^g/ organ ism) of the predator 

k'- = uptake rate of the chemical from the water (L/d-kg(w) 
kg(w) = kg wet wt) 

w = wet weight of organism (kg(w)) 

a- | 1 = chemical assimilation efficiency (fig absorbed//ig ingested) 
1,1 i 

C- ■_■, = specific consumption rate (kg(w) prey/kg(w) predator-d) 
1,1 i 

u'_. = chemical body burden of prey (/ig/organism) , and 

K = chemical excretion rate, d 

Since 

v' = ™ Up (3) 
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where w-,. = organism weight on a lipid basis (fig lipid), Eq, (2) can 
therefore be written using a lipid based concentration of the chemical in 
the organism as follows 

IT - "rf* ♦ «M-1 «!,« "M - 0i ♦ "I* "i (4) 

where now all concentration and parameters (uptake and food consumption 

rates) are expressed on a lipid basis. Also, G = the net growth rate of 

the organism and is given by 

G = aC - r (5) 

where r is the organism respiration rate (d ) and a is the food 

assimilation efficiency. Equation (5) results from a simple growth 

equation on the weight of the organism given by 

| = (aC - r) w = Gw (6) 

The food consumption rate is then estimated from Eq. (5) as 

C = £-£* . (7) 

1. Uptake From Water, Variable Efficiency. - The uptake rate in the 
above equations, k • represents the volumetric rate at which the chemical 
is transferred across the gill membrane of fish, for example. 
Measurements of the efficiency of transfer across the lipo- protein 
membranes of fish have been made (see, for example, (McKim et al. 1985)). 
In general, these laboratory data indicate a reduced transfer efficiency 
at low K values (< 5) , increasing efficiency values at intermediate 
levels of the octanol water partition coefficient and then declining 
values of efficiency at higher K levels of greater than about 7. 

Recognizing that the uptake rate is directly related to the 
respiration rate of the organism and the efficiency of transfer, Thomann 
(1987a) indicates that an approximate relationship for the uptake rate as 
a function of the efficiency of transfer is given by 
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K - 1000 S±£ E (K ou ) (8) 

where E (K ) is the efficiency of transfer as a function of the octanol 

water partition coefficient. This functional relationship is non-linear 

as noted previously and to first approximation for larger organisms of 

weight greater than about 100 g(w) is given by the following equations. 

For log K Qw =2-3: 

log E = -1.5 + 0.4 log K (9a) 

For log K Qw = 3 - 6: 

E = 0.5 (9b) 

For log K 6-10: 

log E = 1.2 - 0.25 log K (9c) 

uw 

These equations represent an empirical relationship between observed 

efficiencies and the octanol water partition coefficient. As such, a full 

mechanistic explanation for the behavior as given by these equations 

requires more detailed study. However, the basic rationale is that at 

lower K values, the transfer across the membrane is rapid through 
ow 

aqueous diffusion layers, but is hindered by the lipid membrane because of 

low fat solubility. At increasing octanol water partition coefficient 

levels, the resistance form the membrane is reduced and the transport is 

proportional to K . At very high K values, water solubility of the 

chemical limits transport and the efficiency decreases. 

The Field BCF 

Using the above expressions for the uptake rate and assuming that at 

equilibrium, the lipid based bioconcentration factor (BCF) is equal to the 

octanol water partition coefficient, the excretion rate can be calculated 

as a function of the K . Under field conditions, however, the growth of 

ow 

the organism must be recognized. That is, in Eq. (4), the relevant "loss 
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mechanism" is the quantity K + G. It can then be shown under field 

generic growth conditions using 

G a 0.01 w -0 - 2 

as representative of the growth rate, that the field BCF (N ) is 

approximately 

N W = K ow I 1 + ElK^f } (U) 

Figure 3.1 from (2) shows that the calculated field log BCF reaches a 

maximum of 5.5 at a log K ., of about 6. At higher log K values, the 

** ow *** ow 

field BCF decreases due to decreased efficiency of transfer and increasing 
effect of organism growth. 

The implication from Fig. 3.1 is that if comparisons are to be made 
between field observed concentration factors and the octanol water 
partition coefficient to determine the degree of bioaccumulation, if any, 
then the comparisons should be made to the estimate field BCF with the 
growth rate of the organism included. 

When this comparison was made for the field observed concentration of 
PCB in the lake trout of Lake Michigan, it was found that the observed 
concentration was about one order of magnitude higher than the 
concentration that would have occurred if only uptake from the water were 
considered (see Thomann, 1987b for further details). 
Phytoplankton BCF 

The above discussion applies to the upper trophic levels of the food 

chain. The laboratory BCF for the base of the food chain, the 

phytoplankton, have been compiled (Thomann 1987b) and are shown in Fig. 

3.2. It appears from these data that the phytoplankton BCF is 

approximately equal to the octanol water partition coefficient. There are 

some indications that the BCF is a constant at log K rt values above about 

ow 
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6. Both relationships are shown in Fig. 3.2. The food chain calculations 

summarized below are particularly sensitive to the BCF of the 

phytoplankton. 

Generic Food Chain Model - Bioaccunulation 

The additional chemical input due to consumption of contaminated prey 
as shown in Eq. (2) can now be included. A steady state generic framework 
is used where the use of "steady state" in food chain compartment models 
has previously been discussed (Thomann 1981). 
Steady State Equations for Field BAF 

The following set of equations can be derived under steady state and 
using a generic four level food chain structure: 

N ul = *jC (12a) 

N 2 . N u2 * f 21 N ul (12b) 

N 3 " N »3 * f 32 N w2 * f 32 f 21 N wl ( 12c) 

N 4 - N w4 * f 43 N w3 * f 43 f 32 N »2 + f 43 f 32 f 32 f 21 N wl < 12d > 

where 

f i,i-i ■ H'.V" 1 (,2e) 

In Eqs. (12), N - represents the field BCF of the ith level of the 
food chain, i.e. 

The cumulative effect of the transfer of chemical from the lower 
levels of the food chain to the top predator level can be seen. Also the 
significance of the phytoplankton BCF, N ■ can be seen. 

The only remaining parameter to be specified is the chemical 
assimilation efficiency, a. This parameter has been measured in the 
laboratory for some chemicals and some organisms. However, additional 
effort to obtain measurements of the chemical food assimilation efficiency 
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is necessary especially over a broad range of octanol partition 
coefficients and food and organism types. 

On the grounds that the mechanism of the assimilation of the chemical 
through ingestion of food is similar to the transfer of the chemical 
across the gill membrane, to first approximation it is assumed that the 
relationship between the assimilation efficiency with log K Qw is given by 
Eqs. (9). 
Comparison to Field Data 

A determination can be made about the adequacy of the model given in 
Eqs. (12) by assigning some generic measures of growth, respiration (and 
resulting food consumption rates) and approximate average weights of each 
trophic level. For simplicity, the fraction lipid and the food 
assimilation efficiency were constant across the food chain at 0.1 and 
0.8, respectively. The average weights for each level of the food chain 
were 0.1, 10, and 1000 g(w) for levels #2, 3, and 4, respectively. 

The BAF for level #4, the top predator is compared to observed BAF 
data in Fig. 3.3. The data are largely from the Great Lakes and in some 
cases required assignment of 7. lipid or the fraction of the water column 
concentration in dissolved form. 

The model calculations indicate that for level #4, the BAF reaches a 
maximum of more than one order of magnitude higher than the calculated 
field BCF with growth included. For a simple no-growth lipid 
partitioning, the field BAF is about one order of magnitude higher at a 
log K of about 6.5. 

To first approximation, then, the generic food chain calculation 
indicates that biomagnif ication is significant only in the region of log 
K between about 5-7. Below log K of about 5, the intake of chemical 
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from water and food is balanced by the loss of chemical due to excretion 
and growth. Above this level however, the excretion rate declines and at 
steady state a net accumulation of chemical occurs which results in 
concentrations above the expected level due to uptake from water only. 
Above log I of about 7, the calculation indicates that the BAF levels 
off and then declines due to decreased efficiency of transfer across the 
membranes of the organism. The model results however, shown in Fig. 3 are 
sensitive to the choice of parameters, specifically the growth rate of 
each trophic level and the BCF of the phytoplankton. 
Sensitivity to Phytoplankton BCF 

The sensitivity of the top level BAF to two estimates of the 
phytoplankton BCF as shown in Fig. 3.2 is shown in Fig. 3.4. As 
indicated, the effect of the constant phytoplankton BCF is to more rapidly 
reduce the BAF at high K Qw to the calculated field BCF. The effect is 
negligible for- log K Qw of < 6. Therefore, if the phytoplankton BCF is 
constant, the BAF for the top level of the food chain at high octanol 
partition coefficients would not ref elect as substantial a food chain 
accumulation as with the phytoplankton BCF equal to the K Qw . 
CONCLUSIONS 

This summary of a food chain model using a generic food chain 

structure indicates that: 

1. For organic chemicals with log K Qw of less than about 5, and for 
chemical screening purposes, complicated food chain modeling is 
not necessary. Rather, a simple lipid partitioning calculation 
using the K as a surrogate for the partition coefficient is 
sufficient. Food chain biomagnif ication effects are apparently 
not significant in this region because of relatively high 
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excretion and growth rates and reduced transfer efficiencies. 

2. For organic chemicals with log K in the range from about 5-7 
food chain magnification may be significant. The degree of 
biomagnification depends on the organism growth rate, the 
phytoplankton BCF and the chemical assimilation efficiency as well 
as other factors. 

The biomagnification in this region is about one order of 
magnitude higher than would be calculated form simple lipid 
partitioning. Therefore, for this group of organic chemicals, 
almost all of the chemical concentration in the upper levels of 
the food chain would be due to consumption of contaminated food 
rather that uptake from the water. 

Also, for this region of log K , a food chain model 
structure is necessary. The degree of complexity of the model may 
vary from a simple steady state generic framework, such as 
discussed in this summary or a more complicated site specific time 
variable, age dependent model. 

3. Above loe K of about 7, model calculations indicate that there 

& ow ' 

is some reason to believe that the field BAF will not continue to 
increase. A plateau and subsequent decline in BAF with increasing 
log K may occur due to decreased transfer efficiency of the 
chemical and growth effects. This region of the model however 
does not have any field data for confirmation and as such the 
model results represent a projection of expected BAF values. 
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FIGURE CAPTIONS 

Fig. 3.1. Comparison of hypothesized field BCF (Eq. 11) for organism 
weight > 10-100 g(w) to simple lipid partition (growth rate = 
0). Uptake efficiency from Eq. (9), p = 0.10. 

Fig. 3.2. Phytoplankton BCF (see [2] for references of data). 

Fig. 3.3. Comparison of calculated BAF as function of K to data (from 
[2], see that reference for data citations). 

Fig. 3.4. Sensitivity of BAF for level #4 to phytoplankton BCF (see [2] 
for references of data) . 
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Figure 3.1 Comparison of hypothesized field BCF (Eq. U) for organism 
weight > 10-100 g(w) to simple lipid partition (growth rate 
- 0). Uptake efficiency from Eq. (9), p - 0.10. 
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Figure 3.2 Phytoplankton BCF (see[2l for references of data) 
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Figure 3.3 Comparison of calculated BAF as function of K to data (from [21, 
see that reference for data citations). 
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ABSTRACT: Correlations between the bioaccumulation tendency of 
organic chemicals in aquatic organisms and 1-octanol-water 
partition coefficients are generally satisfactory for moderately 
hydrophobic chemicals (log K < 5), but show a loss of 

linearity for very hydrophobic compounds (log K > 5) in both 

laboratory experiments and under field conditions. In order to 
establish the possible cause(s) of this phenomenon the roles of 
metabolism, exposure time, membrane permeation, lipid 
solubility, bioavailability, fish growth, chemical uptake from 
food and elimination to the faeces on the bioconcentration 
potential of chemicals and their relative importance are 
discussed. The role of food chain accumulation or 
biomagnif ication is discussed and it is concluded that at each 
trophic level of the food chain a chemical's concentration or 
fugacity increases by a factor of 3 to 5. 

INTRODUCTION 

The ability of organic chemicals to bioaccumulate and elicit a certain 
toxic response in aquatic organisms is usually ascribed to the 
partitioning tendency of these chemicals between water and the organism or 
a particular target site within the organism. Since partitioning of 
organic chemicals is predominantly into the lipid tissue of the organism, 
and 1-octanol resembles these lipids in many respects, various types of 
activities, toxicities and bioaccumulation phenomena have been 
successfully correlated with the 1-octanol -water partition coefficient, 
K . For example linear correlations with respect to K have been 
observed for narcotic action of drugs (Mullins 1954), for acute toxicities 
in fish and crustaceans (Bobra et al. 1985, Veith et al. 1983, Hermens et 
al. 1985, Abernethy and Mackay 1987) and for bioconcentration factors of 
organic environmental contaminants in fish (Neely et al. 1974, Mackay 
1982, Qpperhuizen et al. 1985). Further, transport parameters such as 
kinetic rate constants for uptake and elimination in fish can be related 
to K (Gobas et al. 1986a, Mackay and Hughes 1984, Gobas and Mackay 
1986). However, when K becomes relatively high, a loss of the linear 
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correlation is usually observed (Dearden 1985). In laboratory experiments 
observed bioconcentration factors of high K QW chemicals (i.e., log K Qw > 
5) are usually lower than predicted from the K Qw (Bruggeman 1984, 
Opperhuizen et al. 1985). However, bioconcentration factors of the same 
chemicals in the field, such as in Lake Ontario (Oliver and Niimi 1985), 
are often higher than predicted from K Qw - This severely limits the 
predictive ability of the presently used linear relationships between the 
BCF and the 1-octanol-water partition coefficient and causes uncertainty 
about the actual mechanisms involved. To the class of chemicals with log 
K exceeding 5 belong many persistent contaminants of environmental 
concern such as polychlorinated biphenyls (PCBs) and polychlorinated 
dibenzo-p-dioxins (PCDDs) , some of which have been observed to be 
extremely toxic to mammals and aquatic organisms (Safe 1984) . 

The objective of this paper is to address and investigate some of the 
factors that singly or in combination may cause the apparent loss of 
linear correlation between bioconcentration factors in fish and K Qw . 
These are exposure time, metabolism, bioavailability, membrane permeation 
and lipid-water partitioning. We also discuss briefly the role of food 
uptake and food chain accumulation (or biomagnification) . First, it is 
useful to review the mathematical background of presently used 
bioconcentration models. 
BIOCONCENTRATION IN FISH 

Bioconcentration is the process by which fish or other aquatic 
organisms absorb chemicals only from the water. Bioconcentration does 
therefore not include chemical uptake from food i.e. biomagnification, 
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which we discuss later. 

Since hydrophobic organic chemicals have a relatively low affinity for 
water compared to the fish lipids, a fish when exposed to contaminated 
water absorbs chemicals from the water until the chemical's fungacities in 
water (f ) and fish (fp) are equal i.e. the ratio of chemical's 
concentrations in the fish Cp, and water, C reaches the fish/water 
partition coefficient, Kn (Mackay 1982). 

As has been discussed in more detail elsewhere the fugacity f can be 
viewed as the chemical's "escaping tendency" from a phase; it is linearly 
related to concentration C through a fungacity capacity Z, i.e. C equals 
Z • f (Mackay and Paterson 1982). When the water and fish have reached 
equilibrium and by definition their fugacities are equal, the fish/water 
partition coefficient or bioconcentration factor Kp is Cp/C or Zp/Z . 
For many organic chemicals a simple two compartment (i.e. fish/water) 
model following first order kinetics gives a kinetic description of the 
bioconcentration process. The expression for fish/water exchange assuming 
no variation of the fish volume with time (e.g. no growth) is then given 

V F • Z F • dfp/dt = V L • Z L ■ df L /dt = D F ■ (f w - f L ) (1) 

where the subscripts F, L and W refer to respectively the fish, the fish 
lipids and water, V is the volume (m ) of the respective phases, t is time 
(h) and Dp is the fish/water transport parameter (mol/Pa.h), which is 
related to the more conventional uptake and elimination rate constants k. 
and k« according to: 

k x = D F /(Vp • Z y ) (2) 

k 2 = D f /(V F ■ Z F ) (3) 

Integration of Eq. (1) with a constant fy and an initial fp of zero 
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gives 

f F * f ¥ • (1 -exp(-D F ■ t/V p • Z F )) (4) 

Since F f is C F /Z p and f y is C w /Z y Eq. (5) can be rewritten as: 

C F = C V ' T " (1 ' exp ( " D F ' l / V F ' Z F )} (5) 

V 

Equation (5) is equivalent to the solution of the more traditional 
differential equation describing uptake and release in fish i.e. 

W - k l • C V " k 2 ' C F <6) 

which after integration gives: 

c f = c w • r • (1 " exp (_k 2 ' t} > (7) 

The bioconcentration factor, K c , is defined as the ratio of chemical 
concentration in the fish and water at steady state (i.e. t = ») and thus 
according to Eqs. (5) and (7): 

K c = C F /C W = k 1 /k 2 = Z p /Z v = L • Z L /Z V = L ■ C L /C y = L • K L (8) 
where L is the volumetric lipid content of the fish i.e. V i/Vp- 
CORRELATIONS BETWEEN K Qw AND K c 

Since the lipid tissue of the fish is the principal site for 
bioaccumulation and 1-octanol is often a satisfactory surrogate for 
lipids, linear correlations are usually observed between K^ and K Qw . The 
partition coefficients K Qw and K c can be expressed in thermodynamic 
quantities as the ratio of molar volumes, v, of the respective phases and 
the chemical's activity coefficient y, defined on a Raoult's basis in the 
two phases: 

K = -H Jt (9) 

ow v Q ■ y Q 

K - L • Vy ' YW (10) 

The activity coefficient y expresses the chemical's affinity for a phase. 
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It is inversely related to the chemical's mole fraction solubility x 
(mol/mol) in that phase i.e. y is 1/x for liquid chemicals or y is F/x for 
solid chemicals, where F is the fugacity ratio or ratio of solid to liquid 
solubilities and x/F is the chemical's subcooled liquid solubility (Mackay 
and Shiu 1984). 

From Eqs. (9) and (10) it follows that correlations between K^ and K Qw 
will be linear i.e. WK is constant as long as y^ly Q remains constant. 

For many chemicals, generally with log K Qw between 2 and 5, 
correlations between log K^ and log K 0w are usually linear. However, 
recent experimental data obtained by the authors and others (Opperhuizen 
et al. 1985 and 1986, Burggeman et al. 1984, Muir et al. 1985 and 1986, 
Gobas et al. 1986b) have clearly shown a loss of linear correlation for a 
number of compounds. This is illustrated in Fig. 4.1, where 
bioconcentration factors as listed in Table 4.1 of a series of compounds 
in guppies (Poecilia Reticulata) are plotted versus their experimentally 
determined K . The following conclusions can be reached from Fig. 4.1: 

(1) For compounds with log K s between 2 and 5.5, log K^ and log K Qw 
show generally a linear correlation such as for the data in Table 4.1 
(confidence intervals have a 957. probability): 

log K L = 1.01 (+/" 0.32) ■ log k Qw + 0.323 (+/- 1.5) 

n = 14, r 2 = 0.76 (11) 

(2) Compounds with log K higher than 6 have lower bioconcentration 
factors than expected based on Eq. (11). For these high K Qw compounds 
bioconcentration factors no longer increase with respect to K Qw , but show 
a maximum of log Kr of approximately 6 for compounds with log K s of 5 
then a gradual decrease with a further increase in K . 

(3) PCDDs have much lower BCFs than expected from their high K Qw s (Muir 
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1986, Opperhuizen et al. 1986, Gobas et al. 1986b). It is remarkable, 

however, that this does not seem to be the case for the congener 2, 3, 7, 

8 TCDD, which has been observed to bioconcentrate in rainbow trout (Kq = 

9270) (Branson et al. 1986). 

(4) Other chlorinated and brominated compounds (Opperhuizen et al. 1985, 

Bruggeman et al. 1984, Zitko and Hutzinger 1976) and also some dyes 

(Anliker and Moser 1986) have been observed to have an extremely low or 

perhaps negligible bioconcentration potential. 

EXPOSURE TIME 

It has been shown by several authors that the time needed to reach 

equilibrium between the water and the fish increases with increasing K Qw . 

This is mainly the result of the longer half-life times t^ i.e. In2/k 2 

for compounds with higher K s due to an increased resistance to chemical 

transport in aqueous phases of the fish (Gobas et al. 1986a, Gobas and 

Mackay 1986). It is thus expected that super-hydrophobic compounds (log 

K > 5 or 6) , will require a very long time to reach equilibrium. BCFs 

measured in relatively short experiments (10-30 days) are found to be 

lower than they would be for prolonged experiments. As has been pointed 

out by Hawker and Connell (1985) a curvi-linear relationship between K^ 

and K should therefore be expected. The experimental bioconcentration 
ow 

data (Fig. 4.1), however, clearly suggest a "parabolic" behaviour of K L 
with respect to K in guppies with a maximum bioconcentration factor for 
compounds with log K s between 5 and 6. The exposure time effect should 
cause a leveling of the Kr - K slope and not a maximum followed by a 
drop. It thus seems unlikely that a relatively short exposure time alone 
can explain the observed bioconcentration behaviour of the high K Qw 
compounds. 
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The dependence of %*** on K can be estimated and thus taken into 
account by conducting experiments with compounds of variable K Qw and 
fitting the data to an expression of the type 

where A and B are constants (Mackay and Hughes 1984) . The parameters A 

and B are fish-specific and are related to the fish si2e such that smaller 

fishes have shorter half times (Gobas and Mackay 1986). Conducting 

experiments with small fish such as guppies (average weight 0.2-0.5 g) 

thus eliminates the primary cause of the low bioconcentration due to short 

exposure times. 

METABOLIC TRANSFORMATION 

Metabolism clearly reduces the BCFs of compounds by enhancing the 

depuration rates with respect to the uptake rate. This can be shown by 

extending the fish-water two compartment model to include metabolic 

transformation, represented by the transformation parameter D^ and the 

first order metabolic rate constants kp: 

V p • df F /dt = D F • (f y - f p ) - D R • f p (13) 

Upon integration this equation becomes 

Z P D r Dp + D R 

C F - V ZJ ' Df7Tj[ ' (1 - «P < " jpT* * '» ( 14 ) 

Equation (14) shows that at steady state or infinite time the BCF i.e. the 

ratio of chemical concentrations in fish and water is no longer controlled 

solely by the thermodynamic quantities Zp and Zy but also by the relative 

rates of uptake, elimination and metabolic transformation. For example, 

it can be observed from Table 4.2 that some compounds have 

bioconcentration factors K P „ which are lower than expected from their 

u,exp r 

high K s according to equilibrium relationship that Kp equals 
to 0.048 • K (Mackay 1982). This is due to an increase in the total 
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depuration constant (kn + k 2 ) over the expected elimination constant k 2 

(i.e. k 1 ■ K r ), which can be fully accounted for by the actually 
l L,exp 

measured biotransformation rate constant k R . The correlated 
biotransformation rate constant, k^ CQr (i.e. kj • (1/K^, - */^C,cor^ 
is in agreement with the observed biotransformation rate constant, 
suggesting that indeed biotransformation is the cause for the low 
biotransformation factors of these compounds. 

Most of the information that exists about the bioaccumulation 
behaviour of super-hydrophobic compounds comes from bioaccumulation 
experiments performed with poly-halogenated aromatic compounds such as 
PCBs, PBBs, PCDDs and mirex. Present experience of transformation rates 
of organic chemicals in organisms, and in the environment in general, 
indicates that halogenated aromatics have extremely low biotransformation 
rates. However, it is possible that enzymatic systems exist in organisms 
which can break certain characteristic groups in molecules. 

Although we do not rule out the possible role of metabolic 
transformation it seems unlikely that metabolic transformation is the sole 
cause of the low BCFs of some of the extensively halogenated aromatics. 
BIOAVAILABILITY 

Because of sorption to dissolved or particulate organic matter in the 
water phase, the total water concentration of many organic chemicals may 
not always reflect the fraction of chemical that is truly dissolved and 
thus believed to be available for uptake in fish. In such cases the 
concentration in the fish reflects only the available fraction of chemical 
in the water, which may be lower than the total concentration of chemical 
present in the water phase. Several authors have shown that the presence 
of humic acids and other organic matter in the water phase can reduce the 
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BCF significantly (McCarthy 1983, Landrum et al. 1985). This effect is 
likely to become more important when compounds become more hydrophobic, 
since a larger fraction of the chemical can be expected to exist in the 
sorbed state. Reduced bioavailability is therefore likely to contribute 
to the loss of linear correlation between the BCF and K Qw for the high K Qw 
compounds. However, the observation that chemicals such as 
hexabromobenzene (Zitko and Hutzinger 1976), some higher chlorinated 
naphthalenes (Opperhuizen et al. 1985) and PCDDs (Opperhuizen et al. 1985, 
Gobas et al. 1986a, Muir et al. 1986) have low BCFs compared to the more 
hydrophobic (higher K Qw ) octa- and decachlorobiphenyl cannot be readily 
explained by reduced bioavailability. 

Most sorption partition coefficients are satisfactorily correlated 
with K but relatively few data exist for highly hydrophobic chemicals. 
It is possible that the conventional sorption correlations break down for 
such chemicals and that other factors such as molecular shape play an 
important role in determining sorption and bioavailability. At present 
knowledge about the extent to which hydrophobic chemicals are bioavailable 
to organisms in waters with a natural organic carbon level is scarce, but 
is of major importance for making reliable predictions about the 
bioaccumulation behaviour or toxicity of these chemicals. 
MEMBRANE PERMEATION 

It has been suggested that relatively large molecules, such as many of 
the high K compounds, may experience a hindered membrane passage (Tulp 
and Hutzinger 1978, Muir et al. 1985, Kenaga and Goring 1980). In such 
cases uptake and elimination rate constants will be lower than expected. 
However, it is suspected that BCFs of conservative chemicals will not be 
affected by reduced membrane permeation provided exposure times are long 
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enough, since BCFs are inherently independent of kinetic factors. 

It has also been shown (Hackay and Hughes 1984, Gobas et al. 1986a, 
Gobas and Mackay 1986) that most of the resistance for uptake and 
elimination of super hydrophobic chemicals lies in the aqueous phase of 
the fish and not in organic phases such as membranes. This implies that 
the membrane resistance has to increase by many orders of magnitude in 
order to affect uptake and elimination rates. 

Based on the observed lack of bioconcentration of hexabromobenzene, 
hepta- and octachloronaphthalene and octachlorodioxin it has been 
suggested that molecules with certain spatial dimensions are not able to 
cross the membrane-water interface and enter the membrane (Opperhuizen et 
al. 1984 and 1985). These molecules, which were suggested to have minimal 
internal cross sections (MICS) of 0.96 nm or larger, are thus believed to 
have a low or zero bioconcentration factor. From a mechanistic point of 
view it seems conceivable that a cut-off diameter exists for passive 
membrane diffusion. However, what this cut-off diameter exactly is 
remains unclear since in a recent unpublished bioconcentration experiment 
brominated biphenyls with the same MICS as hexabromobenzene were 
accumulated in guppies (Gobas et al., 1987). The obvious test of this 
effect is to perform bioconcentration experiments with conservative 
chemicals varying in their spatial dimensions. 
LIPID/WATER PARTITIONING 

As pointed out earlier, the linear correlation between K^ and K is 
the result of the similarity in phase properties of the fish lipids and 
1-octanol towards organic chemicals i.e. yy /y is constant. It can thus 
be argued that the loss of linear correlation between Kq and K Qw for high 
K compounds may be due to the fact that yT/y is no longer constant. 
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Two arguments support this contention: 

a) the fish lipid phase consists primarily of biological membranes, 
in which the molecules are predominantly arranged in bilayers. The lipid 
phase thus has a distinct structure and restricted spatial dimensions. 
Since the 1-octanol phase is a bulk phase presumably with little or no 
structure, organic solutes may display different activity coefficients and 
partitioning behaviour in 1-octanol than in membranes. 

b) it can be shown that the subcooled liquid solubility in 1-octanol 
i.e. 1/y decreases with an increase in the molar volume of the solute. 
This has been suggested previously (Miller et al. 1985) and is illustrated 
in Fig. 4.2, where x /F (i.e. l/y Q ) of a series of high K Qw compounds is 
plotted versus the molar volume of the solute, calculated according to the 
Le Bas method (Reid et al. 1977). The values for x Q were obtained by 
measuring the solubilities of these chemicals in 1-octanol by a method 
described earlier (Miller et al. 1985) and from Anliker and Moser (1986). 
Since all chemicals were solids at the temperature of interest i.e. 25 C 
the fugacity ratio F was calculated as 

F = exp (AS(1 - T M /T)/R) (15) 

where AS is the entropy of fusion at the melting point and estimated to be 
13.5 calorie entropy units (Yalkowsky 1979), T M (K) is the compound's 
melting point, T(K) is the experimental temperature and R is the gas 
constant i.e. 1.98 cal/mol.K. The experimentally obtained solubilities in 
1-octanol and calculated fugacity ratios are listed in Table 4.3. It 
follows from Fig. 4.2 that all chemicals have mole fraction solubilities 
below 1 and accordingly activity coefficients above 1. This implies that 
these chemicals do not form ideal solutions in 1-octanol and that 
1-octanol has a limited capacity to accommodate these chemicals. Figure 
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4.2 also shows that the log solubility in 1-octanol is approximately 

inversely related to the chemical's molar volume according to 

log y Q = 0.0051 (+/- 0.0008) • V - 0.26 (+/- 0.32) 

n = 41, r 2 = 0.80 (16) 

showing that 1-octanol has a lower capacity to accommodate larger volume 

molecules. 

As has been demonstrated convincingly by Kamlet et al., (1983) and 
Kamlet (1986), the best correlation of solubilities in solvents in general 
is achieved when both molar volume parameters and parameters expressing 
the interactions between solute and solvent are introduced. However, for 
solutes, which in this case are closely related chemically, variations in 
solubilities due to different solute-solvent interaction can be expected 
to be relatively small compared to the influence of the solute's varying 
molar volume. It may thus be concluded that the partitioning tendency of 
organic solutes in 1-octanol is not only determined by the behaviour of 
the solute in the water phase i.e. y^ but also by the behaviour in the 
organic phase. It also seems that for the class of aromatic compounds the 
formation of a cavity in 1-octanol of the dimensions of the solute 
molecule is the most important process determining the solute's subcooled 
liquid solubility. It can thus be expected that this cavity formation 
also governs the solubilities of these molecules in membranes. However, 
when solute molecules become relatively large i.e. have dimensions 
comparable to the thickness of the membrane phase, the structured membrane 
phase may not appear as a homogeneous phase to the solute, and specific 
interactions such as between the polar head groups of the phospholipids 
may have to be broken in order to dissolve the solute molecule in the 
hydrophobic interior of the membrane. It thus seems conceivable that for 
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relatively large solutes cavity formation in the membrane phase is 
relatively more difficult than for smaller molecules. Accordingly the 
subcooled liquid solubilities per unit of molar volume i.e. x/Fv of 
relatively large molecules can be expected to smaller and the activity 
coefficient per unit of molar volume i.e. y/v correspondingly higher than 
for smaller molecules. The expected behaviour of y T and also of y with 
respect to molar volume for solutes, for which membrane/solute 
interactions are approximately equal is illustrated in Fig. 4.3. This 
figure shows that the relationship between K rt and K r cannot be linear but 

OW \j 

must be "parabolic" with a maximum value for Kp. 

Preliminary work in our laboratory on membrane/water partitioning of 
hydrophobic compounds of varying molar volume supports this view. 
Partition coefficients between dymiristoyl-phosphatidylcholine (DMPC) 
membrane vessicles and water of chlorobenzenes, PCBs and PCDDs show a 
linear correlation with respect to the solute's molar volume up to a molar 
volume of 200 cm /mol. A maximum membrane-water partition coefficient of 
10 was observed for solutes with a molar volume of approximately 300 
cm /mol. For larger solutes the membrane-water partition coefficient 
tends to fall. If, as seems likely, artificial membranes do indeed 
simulate lipids in organisms better than 1-octanol then it is likely that 
the loss of linear correlation between Kp and K is a least partly due to 
differences in solvent characteristics between natural lipids and 
1-octanol. It is conceivable that there may emerge a "rule" that for 
molecules less than a certain volume or with certain dimensions, 1-octanol 
is satisfactorily surrogate, i.e. y Q and y L are approximately equal. But 
for larger molecules this similarity breaks down and y, is much larger 
than y Q and 1-octanol is no longer a satisfactory surrogate, and K a 
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satisfactory descriptor. In that case membrane-water partition 
coefficients may be used as a more reliable but still simple method to 
estimate and correlate the bioconcentration factors of organic chemicals 
in aquatic organisms. 
FISH GROVTH 

If the fish is growing the term Vp in Eq. (1) is no longer constant 
and its time dependence must be included. The nature of the effect is 
that growth causes dilution, thus there is an apparent (but unreal) 
increase in k„ or Dy. This effect can be included easily by adjusting Eq. 
(1). Elsewhere in these proceedings Thomann has discussed this issue in 
more detail. 
CHEMICAL UPTAKE FROM FOOD 

Chemical uptake from food and chemical excretion to the faeces can 
also be included in the differential equation describing chemical uptake 
in fish. We do not report here a full account of this issue but we report 
on some recent findings. 

For substances of relatively low K (i-e., log K < 5), the chemical 
exchange between fish and water (via the gills) is relatively fast 
compared to the chemical uptake from food and excretion to the faeces. As 
a result the chemical ingested with the food is rapidly equilibrated with 
the water through gill exchange and the concentration in the fish will 
approach or achieve equilibrium with the water. Therefore it does not 
matter if the fish obtains the chemical from the water or food. In that 
case the bioaccumulation factor due to combined chemical uptake from water 
and food corresponds to the bioconcentration factor measured in laboratory 
experiments and, when not metabolized, can be correlated with the 
1-octanol-water partition coefficient. For these chemicals there is 
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little or no biomagnif ication. 

As K increases, the time needed to release the chemical to the water 
ow ' 

increases. The fish may thus achieve a state in which the rate of intake 
from food exceeds the release rate, which can result in a fish 
concentration or fugacity exceeding the equilibrium value. This 
phenomenon, which is discussed in more detail elsewhere (Gobas et al. 
1987), is generally referred to as biomagnif ication, Ve believe that the 
mechanism underlying this process is that as a result of food digestion in 
the gastro-intestinal tract the organic mass in the Gl-tract continuously 
decreases thus increasing the chemical's concentration or fugacity in the 
Gl-tract and driving the chemical across the gut wall into the fish. 

A fish may then live in water contaminated with a chemical, which has 
a fugacity in the water of for example 10 Pa and eat small organisms 
with a similar chemical fugacity, but achieve a chemical fugacity of 
perhaps 4 • 10 Pa. This increase, in the chemical's fugacity takes place 
at every trophic level of the food chain, i.e., where one organism feeds 
on another organism of a lower trophic level. So, if the fish with a 
chemical fugacity of 4 * 10 is eaten by another organism of a higher 
trophic level, that organism may reach a chemical fugacity of 16 • 10 
Pa. It appears from our work that fish eating birds such as herring gulls 
achieve a chemical fugacity increase of approximately 10 to 30. This 
biomagnif ication phenomenon appears to be of the order of a factor of 3 to 
5 per trophic level for aquatic organisms and predominantly applies to 
chemicals with a log K Qw between 5 and 7. At higher K Qw values we suspect 
that the efficiency of gut absorption efficiency decreases (Gobas et al. 
1987). These chemicals may therefore not biomagnify or even 
bioconcentrate appreciably. 
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It is worth noting that from a regulatory standpoint that this 
phenomenon is desirable, but it is essential to understand its cause since 
it is conceivable that chemicals may exist of high K Qw which do 
bioaccumulate. In addition it should be noted that bioconcentration data 
for chemicals with log K exceeding 7 is extremely scarce. Conclusions 
about their bioaccumulation potential are therefore difficult to draw. 

We are presently working to assemble and validate models describing 
these processes. We focus our attention on chemicals of very high K 
(i.e., log K between 5 and 8) since they are the chemicals with the 
highest expected bioaccumulation potential. 
CONCLUSIONS 

In this paper we have discussed the factors that may affect the linear 
correlation between bioconcentration factors in fish and the 
1-octanol-water partition coefficient. Since a lipid-water partitioning 
process is the principal process controlling the uptake of organic 
chemicals from the water the same factors may also affect other QSARs in 
aquatic toxicology. 

For compounds with log K between 1 and 5 correlations between K Qw 
and Kp are generally linear. When deviations from linear behaviour occur, 
metabolic transformation has to be considered the principal reason for 
this effect, since all other factors influencing K Qw - Kq relationships 
apply only to high K compounds. 

For compounds with log K exceeding 5 or 6 simple linear 

L/W 

relationships between Kp and K breakdown. In laboratory experiments 

involving chemical exchange between water and fish K P often shows a 

maximum for compounds with log K of 5 to 6 and then tends to fall for 

compounds with log K exceeding 6. For these chemicals an equilibrium 
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between fish and water is reached only after a very long exposure time, 
which often exceeds the experimental exposure time, such that it is often 
difficult to determine the actual steady state level in the fish and thus 
the true bioconcentration factor. Indeed the required exposure time may 
exceed the lifetime of the fish. However, even when curve fitting or 
extrapolation procedures are used to determine steady state concentrations 
in the fish, the calculated bioconcentration factors are still lower than 
would be expected from their high K . The combination of a relatively 
short exposure time, slow uptake and the "diluting" effect of fish growth 
does therefore not explain the relatively low bioconcentration factors of 
the high K compounds. 

The decrease in subcooled liquid solubility of organic solutes in 
1-octanol with increasing molar volume shows that the partitioning 
tendency of these chemicals is not only determined by their behaviour in 
the water phase i.e. their hydrophobicity or y,, but also by their 
behaviour in the organic medium. It is conceivable that differences in 
phase characteristics between 1-octanol and the organism lipids towards 
organic solutes affect the relationship between 1-octanol-water and fish 
lipid/water partition coefficients. Preliminary studies on membrane/water 
partition coefficients indicate that lipid phases in organisms, which 
predominantly consist of membranes, may have a relatively low potential to 
accommodate larger molar volume molecules, many of which have high K . 
It seems therefore conceivable that the relatively low solubility of 
voluminous molecules in membranes compared to 1-octanol is one of the 
causes of the loss of linear correlation between Kp and K . 

A second reason for the loss of this linear correlation can be a 
reduced bioavailability of high K compounds. Unfortunately experimental 
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data which can distinguish between the bioavailable and non-bioava liable 
fraction of (super) hydrophobic chemicals in water are scarce because of 
the extreme experimental difficulty of making these measurements. As a 
result, no accurate predictions can yet be made about the role of a 
reduced bioavailability on the linearity of correlations between K^ and 

K ow' 

Under field conditions, where a fish is exchanging chemical between 

water and food, correlation between the bioaccumulation factor and the 

1-octanol-water partition coefficient also show a loss of linear 

correlation. For high IT chemicals bioaccumulation factors are often 

ow 

higher than would be predicted from the 1-octanol-water partition 
coefficient. This is probably the result of biomagnif ication which can 
result in fish-water concentration ratios which exceed the fish-water 
partition coefficients determined in the laboratory. 

We conclude that when assessing the behaviour of new and existing 
chemicals it is unwise to extrapolate existing K Qw - K^ correlations into 
the super-hydrophobic region. It is believed that a chemical by chemical 
assessment will be needed in which the roles of metabolism, 
bioavailability, membrane permeation, lipid/water partitioning, fish 
growth and food uptake are considered as modifying the basic equilibrium 
relationship. Some of these phenomena are fortunate in that they mitigate 
the bioaccumulation tendencies of many chemicals but, since these 
chemicals have a very high bioaccumulation potential and are potentially 
very toxic, it seems unwise to exploit or depend on these factors without 
fully understanding their nature. The phenomenon of food chain 
accumulation has to be examined carefully since it is likely to result in 
very high concentrations in organisms of higher trophic level. 
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Finally, we believe that progress towards assembling a comprehensive 
model of chemical uptake in fish, including equilibria, kinetics and food 
chain effects (i.e., biomagnification) is best achieved by assembling 
models of various types, testing them and modifying them to improve the 
fidelity with which they describe bioaccumulation. 
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TABLE 4.1: Chemical Abstract Nuabers (CAS), Octanol /Water Partition Coefficients 
(log kow) BCFs In Whole Fish (log kc) and in lipids (log Kl) in the 
guppy (Poecilia reticulata} 

COMPOUND 

2,5 DLCHLOROBIPHENYL 

2,2,5,5 TETRACHLOROBIPHENYL 

2, 2, 4, 1,5, 5 HEXACHLOROBIPHENYL 

2,2,3,3,1,4,5,5, OCTACHLOROB I PHENYL 

DECACHLOROBIPHNEYL 

HEXABROMOBENZENE 

2 -MO NOCHLORO NAPHTHALENE 

1,4 OICHLORONAPHTHALENE 

1,8 OICHLORONAPHTHALENE 

2,3 OICHLORONAPHTHALENE 

2,7 OICHLORONAPHTHALENE 

1,3,7 TRICHLORONAPHTHALENE 

1,2,3,4 TETRACHLORONAPHTHALENE 

1.3.5.7 TETRACHLORONAPHTHALENE 

1.3.5.8 TETRACHLORONAPHTHALENE 



-j PENTACHLOROBENZENE 

m 



2,3',4',5 TETRACHLOROBIPHENYL 

1,4 0ICHL0ROBENZENE 

1.2.3 TRICHLOROBENZENE 
1,3,5 TRICHLOROBENZENE 
1,2,3,5 TETRACHLOROBENZENE 
PENTACHLOROBENZENE 
HEXACHLOROBENZENE 

2,7, DICHLOROOIBENZO-P-DIOXIN 

1.2.4 TRICHLORODIBENZO-P-DIOXIN 
1,2,3,4 TETRACHLORODIBENZO-P-DIOXIN 
OCTACHLORODIBENZO-P-DIOXIN 

a: froa Shiu and Hckay 

b: froa Opperhulzen et al . 19S5 

c: fro* Miller a_t a±. 1985 

d: froa Bruggeraan e t a 1 ■ 1984 

e: froa Koneaann et a 1 , 1980 

f: froa Gobas £t^ _aj_. l9B6b 

g; froa Andren st_ a_K 1987 



CAS Nuaber 


Log Kow 


Log Kc 


Log Kl 


34883-39-1 


5.10 


a 


4.26 


5.72 d 


35693-99-3 


6.10 


a 


4.92 


6.38 d 


35065-27-1 


6.90 


a 


5.32 


6.78 d 


35694-08-7 


7.10 


a 


4.35 


5.81 d 


2051-24-3 


8.26 


a 


4.02 


5.48 d 


87-82-1 


> 6.0 




< 5.0 b 


91-58-7 


3.90 


b 


3.63 


4.81 b 


1825-31-6 


4.66 


b 


3.36 


4.54 b 


2050-74-0 


4.19 


b 


3.79 


4.96 b 


2050-75-1 


4.51 


b 


4.04 


5.22 b 


2198-77-8 


4.56 


b 


4.04 


5.22 b 


1321-65-9 


5.35 


b 


4.43 


5.61 b 


20020-02-4 


5.75 


b 


4.66 


5.71 b 


53555-64-9 


6.19 


b 


5.06 


5.71 b 


31604-28-1 


5.81 


b 


4.69 


5.57 b 


608-93-5 


5.03 


c 


4.38 


5.40 b 


70632-48-0 


5.90 


a 


4.92 


5.50 b 


106-46-7 


3.38 


c 


1 .99 


3.25 • 


87-61-6 


4 .04 


c 


2.84 


4.11 e 


108-70-3 


4.02 


c 


2 .88 


4.14 • 


634-90-2 


4.65 


c 


3.59 


4.85 a 


608-93-5 


5.03 


c 


4.14 


5.41 a 


118-74-1 


5.45 


c 


4.19 


5.46 a 


33857-26-0 


5.75 


9 


1 .7 


2.9 f 


39227-58-2 


6.45 


g 


1.9 


3.1 f 


30746-58-8 


6.70 


9 


2.2 


3.4 f 


3268-87-9 


8. 40 


9 


1 .9 


3.1 f 



Table 4.2: Experimentally observed (exp) and correlated (cor) bioconcentration 
factors (Kc) and metabolic rate constants (k R ) (1/hr) in fish. 



COMPOUND log Kow Kc.tip Kc,c«r ki ki,«i P k«,cor 



ANTHRACENE { 1 ) 
BENZO-a-PYRENE (1) 
ACRIDINE (2) 
BENZ[a]ACRIDINE (2) 
DIBENZ[a,h]ACRIDINE (2) 



4.4 5 


675 


1350 


36 


0.027 


0.027 


6.06 


490 


55100 


49 


0.31 


0.1 


3.3 


100 


96 










4.45 


106 


1350 


100 


0.49 


0.87 


5.6 


107 


19100 


12 


0.18 


0.11 



(1) Spacie et al . , 1983 

(2) Southworth et al . , 1980 
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TABLE 4.3: CAS numbers, molecular weight, aelting point (Te), Molecular weight (Hw), fugaclty ratio (F), solubility in 
1-octano (Cs), activity coefficient (log yo), and 1-octanol/water partition coefficient (log kow) 

COMPOUND 

4 Monochlorobiphenyl 
4.4 dichlo rob I phenyl 
2,4,5 trichlorobi phenyl 
2,2,5,5 tetrachlorobiphenyl 

2.3.4.5 tetrachlorobiphenyl 
decacn 1 or ob I pheny 1 
naphthalene 
bi phenyl 
scenaphthene 
fluorene 

1 -aethy I -f I uorene 
anthracene 
phenathrene 
f luoranthene 
pyrene 
chrysene 

2.3 benzofluorene 
perylene 

1.2.5.6 d i banzanthracene 

Fro* Mosar and An li leer: 

DDT 

hexach 1 o rocy c 1 ohexane 

hexach I orobenzene 

1.4 dfchlorobenzene 
1,2,4 trlchlorobenzene 

CIBA-GEICY LIPOPHILIC DYE #: 

I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

XII 

XIII 

XIV 

XV 

XVIII 

XIX 

XX 

STRUCTURAL FORMULAS FOR THE LISTED LIPOPHILIC DYES CAM BE FOUMD IN ANLIKER AND HOSAR (1986). 



CAS Nuaber 


MOL.W. 


TM (K) 


V(ca 


3/aol) 


F 


C(s)m1/1 


Log yo 


Log Ko* 


2051-62-9 


188.7 


351 




206 




2.08E-01 


6.08E-01 


0.56 


4.50 


2050-68-2 


223.1 


422 




226 




5.46E-02 


7.13E-02 


0.78 


5.30 


15862-07-4 


257.5 


350 




247 




2.92E-01 


1.76E-01 


1.12 


5.60 


35693-99-3 


292 


360 




268 




2.28E-01 


2.34E-01 


0.69 


6.10 


33284-53-6 


292 


364 




268 




2.09E-01 


1.40E-01 


1.07 


5.90 


2051-24-3 


467.7 


579 




394 




1.48E-03 


1.69E-03 


0.84 


8.26 


91-20-3 


128.2 


353 




148 




2.68E-01 


9.59E-01 


0.34 


3.35 


92-52-4 


154.21 


343 




185 




3.42E-01 


6.91E-01 


0.59 


3.95 


83-32-9 


154.21 


369 




173 




1.B5E-01 


2.56E-01 


0.76 


3.92 


86-73-7 


154.21 


369 




188 




1.17E-01 


2.24E-01 


0.62 


4.18 


1730-37-6 


166.21 


358 




210 




2.39E-01 


2.96E-01 


0.80 


4.97 


120-12-7 


178.2 


469 




197 




1.17E-02 


1.18E-02 


0.89 


4.54 


85-01-6 


178.2 


374 




199 




1.66E-01 


4.01E-01 


0.51 


4.57 


206-44-0 


202.3 


384 




217 




1.32E-01 


1.73E-01 


0.78 


5.22 


129-00-0 


202.3 


429 




214 




4.67E-02 


1.40E-01 


0.42 


5.18 


218-01-9 


226.3 


528 




251 




4.78E-03 


2.00E-03 


1.26 


5.79 


243-17-4 


216.28 


482 




240 




1.38E-02 


1.77E-02 


0.79 


5.75 


196-55-0 


252.3 


550 




263 




2.8SE-03 


3.04E-03 


0.67 


6.50 


53-70-3 


278.35 


539 




300 




3.71E-03 


9.38E-04 


1.49 


7.19 


50-29-3 


354.48 


382 




364 




1.37E-01 


8.18E-02 


1.12 


6.36 


319-84-6 


290.82 


433 




259 




4.24E-02 


5-16E-02 


0.81 


3.80 


118-74-1 


264.76 


503 




221 




8.46E-03 


1.53E-02 


0.64 


5.47 


106-46-7 


147 


326 




138 




4.99E-01 


1.77E*00 


0.35 


3.38 


120-82-1 


181.45 


290 




159 




1 .00E+00 




0.00 


3.98 




546.38 


448 




565 




3.00E-02 


1 .48E-04 


3.20 


2.50 




333 


496 




349 




9.49E-03 


3.00E-04 


2.40 


3.40 




363 


498 




422 




9.49E-03 


3.31E-04 


2.36 


>4.0 




366.19 


492. 


2 


364. 


6 


1.09E-O2 


4.62E-04 


2.27 


>4.0 




446.53 


421. 


2 


552. 


3 


5.59E-02 


9.14E-04 


2.68 


4 




423.3 


457. 


2 


541. 


6 


2.44E-02 


1.69E.03 


2.01 


>2.0 




426.9 


443 


2 


459. 


5 


3.37E-02 


2.23E-03 


2.08 


« 




528 


390 


2 


528. 


3 


1.14E-01 


3.16E-03 


2.46 


4.1 




373.83 


419 


2 


428. 


7 


5.86E-02 


5.62E-03 


1.92 


3.9 




411.54 


446 


2 


504 




3.14E-02 


5.90E-03 


1.62 


4.5 




340.33 


549 


2 


399. 


2 


2.93E-03 


2.76E-05 


2.92 


3.82 




307.29 


563 


2 


366. 


7 


2.12E-03 


5.53E-05 


2.46 


5.4 




486.76 


563 


2 


S46. 


2 


2.12E-03 


1.60E-05 


3.02 


8.38 




629.5 


593 


2 


735 




1.06E-03 


7.94E-07 


4.02 


6.6 




726.44 


593 


2 


623. 


2 


1.06E-03 


6.33E-07 


4.12 


8.1 




439.76 


603 


2 






8.45E-04 


7.96E-06 


2.92 


3.4 




312.32 


673 


2 


342. 


2 


1.68E.04 


5.44E-06 


2.39 


4.3 




1100 


753 


2 


960 




2.673-05 


6.36E-08 


3.52 


17.4 



11267-108.1 



FIGURE CAPTIONS 

Fig. 4,1. Bioconcentration factors based on fish lipid weights of 
chlorobenzenes (ClB), PCBs, polychlorinated naphtalenes (PCN), 
chlorinated dibenzo-p-dioxins (PCDD) in Poecilia reticulata 
(guppy) as a function of K . 

Fig. 4.2. Activity coefficients of chlorinated aromatics, polynuclear 
aromatics and lipophilic dyes in 1-octanol as a function of 
molar volume. 

Fig. 4.3. The expected behaviour of the activity coefficient of 
chemically related compounds in water (V), 1-octanol (0) and 
lipid tissues of organisms (1) versus molar volume, v. 
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Figure 4.1: 



Bioconcentration factors based on fish lipid weights of chlorobenzenes 
(C1B), PCBs, polychlorinated naphthalenes (PCN), chlorinated 
dibenzo-p-dioxins (PCDD) in Poecilia reticulata (guppy) as a function 
of Kow. 
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Figure 4.2: Activity coefficients of chlorinated aromatics, polynuclear 

aromatics and lipophilic dyes in 1-octanol as a function of molar 
volume. 
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Figure 4.3: The expected behaviour of the activity coefficient of chemically 
related compounds in water (W), 1-octanol (0) and lipid tissues 
of organisms (1) versus molar volume, v. 
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INTRODUCTION 

The benthic food chain is an important component of the transport and 
ecosystem fate of organic chemicals in the Great Lakes. Benthic animals 
can significantly alter the physical and chemical properties of sediments, 
thus affecting constituent mass fluxes across the sediment-water 
interface. Sediment exposures can cause accumulation of organic chemicals 
in benthic animals, with the potential for subsequent trophic transfer to 
higher predators. 

The impetus for this work was the desire of the Ontario Ministry of 
the Environment to incorporate an aquatic food chain component into 
existing transport and fate models for organic chemicals in the Great 
Lakes. The goal of this overall effort is to identify essential state 
variables and process mechanisms that govern the fate of these chemicals 
in aquatic ecosystems. Most of the existing research in this area has 
been on the pelagic food chain. Comparatively little is known about 
organic chemical dynamics in the benthic food chain. 
Scope 

The focus in this paper is on the benthic component of the overall 
aquatic food chain. Emphasis is placed on the effects of sediment 
exposures on benthic animals, and on the potential for trophic transfer to 
higher predators. Results from this work have direct implications for the 
degree of model complexity required to describe organic chemical dynamics 
and fate in the benthic food chain. For completeness, a brief review of 
pertinent literature is presented on the effects of benthic animals on the 
physical-chemical properties of sediments. Results from this review may 
have implications for the physical-chemical component of the overall 
ecosystem model. 
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STRUCTURE OF BENTBIC FOOD CBAIN 

There is considerable uncertainty regarding the structure and function 
of the benthic food chains in the Great Lakes. Flint (1986) has recently 
developed a comprehensive conceptual model of the Lake Ontario deepwater 
food web. The model is based on a synthesis of existing information on 
standing stocks, production rates and potential trophic linkages. Figure 
5.1 contains a schematic diagram of the hypothesized trophic relationships 
and Fig. 5.2 contains a schematic analysis of annual carbon production and 
energy flow in this food web. Oliver and Niimi (1987) have proposed a 
more simplified model of the Lake Ontario food chain (Fig. 5.3). 

The principal differences between these two models are in the coupling 
between the benthic and pelagic food chains. The Flint model contains 
direct linkages between macrobenthos (oligochaetes, chironomids and 
amphipods) and adult alewives, smelt, and large fish (primarly salmonids) 
in the pelagic food chain. The Oliver and Niimi model does not contain 
these linkages. Furthermore, the Flint model links sculpins (benthic 
fish) to smelt (pelagic food chain), whereas the Oliver and Niimi model 
links sculpins directly to large fish. Despite uncertainties in actual 
trophic linkages, it is clear that the main food source for the benthic 
community is from phytoplankton photosynthesis and phytoplankton-derived 
detritus from the overlying waters. 

The structure of the benthic food chain has implications for the types 
of models required to describe the distribution of contaminants. Thomann 
and Connolly (1984) and Oliver and Niimi (1987) have shown that chemical 
concentrations in the pelagic food chain depend on the trophic position of 
the animal, especially the number of trophic transfers in the food uptake 
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pathway. In the benthic food chain simple empirical models may suffice 
for animals at low trophic levels, while mechanistic food chain models may 
he required for animals at higher trophic levels. 
EFFECTS OF BENTHOS ON PHYSICAL-CHEMICAL PROPERTIES OF SEDIMENTS 

Benthic organisms can have significant effects on the physical and 
chemical properties of sediments and hence, on the entire sediment 
exposure environment. Lee and Swartz (1980) have reviewed the mechanisms 
and rates of biodeposition and bioturbation in the marine environment. 
They concluded that the major biogenic processes modifying sediments are 
biodeposition, particle transport, fluid transport, stabilizing mechanisms 
and macrofaunal-microbial interactions. In their opinion, there is 
currently a basic understanding of biodeposition and particle transport, 
but not of the latter three processes. 

McCall and Tevesz (1982) have reviewed the effects of freshwater 
benthos on physical properties of bottom sediments. They focused on 
macrobenthos (adult length > 1 mm) because, in their opinion, too little 
is known about freshwater meio- and microbenthos. McCall and Fisher 
(1980) studied the effects of tubificid oligochaetes on Lake Erie 
sediments and concluded that the oligochaetes could alter grain size, 
settling velocity, erodibility, porosity, permeability, vertical gradients 
of particulate and dissolved sedimentary materials, diffusive fluxes and 
sediment oxygen demand. Observed variations in these parameters ranged 
from factors of 2 to 10 in the presence of the tubificid oligochaetes. 
McCall et al. (1979) studied sediment reworking activities of unionid 
bivalves in Lake Erie and found that they plow the top 10 cm and they pump 
large quantities of water through the sediment. These filter-feeding 
clams are the largest invertebrates in Lake Erie and they account for 
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most of the secondary production in the western basin. 

Rob-bins (1986) developed a model to describe the processes of sediment 
diagenesis and the effects of biogenic reworking by conveyer-belt, 
deposit-feeding benthic organisms (e.g., tubificid oligochaetes) . The 
model is based on transport and chemical reaction and it includes 
compaction, diffusion and advective redistribution by the benthic 
organisms. Robbins concluded that particle diffusion, as well as 
a.dvection, was required to accurately describe the distributions of highly 
sorbed tracer constituents. 

Karickhoff and Morris (1985) used laboratory microcosms to study the 
transport of hydrophobic pollutants in bottom sediments in the presence of 
oligochaete worms, as contrasted to transport by molecular diffusion. 
They found that in the presence of stable worm populations, pollutant 
transport within the biologically mixed surface layer was entirely worm 
mediated. More than 90 percent of the chemicals in this layer were 
transported to the sediment surface in a 30-50 day period. Sediment was 
vertically pumped to the surface at rates from 9 to 30 times faster than 
by maximum molecular diffusion rates for the highly sorbed chemicals used 
in the experiments. Mass flux from sediment to water column, however, was 
only enhanced by a factor of 4 to 6. The mitigating process was sediment 
pelletization which trapped sorbed chemicals and significantly retarded 
chemical release into the aqueous phase at the sediment-water interface. 

Gschwend et al. (1986) have proposed a modeling framework to predict 
the transfer of hydrophobic organic chemicals between bed sediments and 
overlying waters. The principal components in this framework are a 
description of the microscopic scale process of sorption kinetics and 
models of the exposure of bed particles to adjacent waters of varying 
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composition. Effects due to diffusion of solutes in interstitial fluids, 

biological mixing of surficial sediments and suspension of bed solids are 

included. 

BIOACCIMULATION BY AQUATIC ORGANISMS 

Environmental contaminants can be accumulated by fish and benthic 
organisms through direct uptake from water and through food chain uptake. 
The primary mode of accumulation in a particular case depends on the 
properties of the chemical, the relative concentrations of the chemical in 
the water and food, and on the trophic position of the animal. 

The term bioconcentration refers to accumulation of a chemical in fish 
by transport through the gill or other membranes, and it excludes 
accumulation by the consumption and absorption of contaminated food (Veith 
et al. 1979). Bioconcentration is usually expressed in terms of a 
bioconcentration factor (BCF) which is the ratio of the chemical 
concentration in the fish to the chemical concentration in the water under 
equilibrium conditions. The term biomagnif ication refers to accumulation 
of chemicals through contaminated food. 

Bioconcentration is generally viewed as a simple phase partitioning 
process. Neely et al. (1974) observed that the bioconcentration of 
organic chemicals in trout followed a linear relationship with the 
octanol-water partition coefficient (K ) of the chemicals. Veith et al. 
(1979) extended this concept to a larger range of organic chemicals using 
uptake experiments with fathead minnows. They pointed out that 
laboratory-derived BCFs were only a conservative estimate of actual 
environmental residues because these residues are the result of potential 
biomagnification, as well as bioconcentration. 
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Mackay (1982) examined the physical-chemical factors influencing 
bioconcentration of organic chemicals and concluded that there was a sound 
theoretical basis for expecting a good correlation between BCF and K Qw . 
His arguments were based on the assumption that an organism, such as a 
fish, could be viewed as consisting of one or more phases of differing 
chemical composition and volume fraction. At equilibrium the chemical 
potential, or fugacity, of the bioconcentrating chemical will be equal in 
the exposure medium and in each phase. 

Chiou (1985) also demonstrated that K Qw was a good predictor for 
lipid-water partition coefficients and for BCFs in fish. He pointed out 
that discrepancies can occur between BCF and K under the following 
conditions: first, chemicals that are unstable in water or that are 
readily metabolized by animals will not reach an equilibrium state; 
second, BCF values obtained with exposure times that are too short may not 
represent true equilibrium between the biotic and water phases; and third, 
macromolecular organic molecules (e.g., humic substances) could enhance 
the apparent solubility of relatively insoluble chemicals, thus decreasing 
the computed BCF value. 

Several studies have demonstrated the importance of exposure levels 
and times to equilibrium in developing valid correlations between BCF and 

K . Oliver and Niimi (1983) studied bioaccumulation of ten 

ow * ' 

chlorobenzenes by rainbow trout over a 119-day exposure period in the 

laboratory. They observed that time to equilibrium was dependent on the 

exposure level and the degree of chlorine substitution. Good correlations 

were obtained between BCF and K„ for chemicals that reached equilibrium 

ow 

(log K < 5) ; however, hexachlorobenzene only obtained equilibrium at a 

L/VV 

high exposure level, not at a lower exposure. Good agreement was obtained 
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using their laboratory-based BCFs to predict fish concentrations in Lake 
Ontario from ambient concentration levels, except for hexachlorobenzene. 
In a subsequent study, Oliver (1984) found that concentration factors for 
24 organic chemicals in oligochaete worms continued to increase over a 
110-day exposure to naturally contaminated Lake Ontario sediments. 
Inconsistent relationships were observed between the BCF and K values 
for the different chemicals. 

There are considerable reported differences in the relative importance 
of water vs food uptake even for organisms at the base of the benthic food 
chain. For example, Adams et al. (1985) studied kepone uptake by 
chironomids and concluded that bioaccumulation occurs primarily via 
interstitial water and not from ingestion of sediments, food or exposure 
to the water column. Lynch and Johnson (1982) studied availability of 
radiolabeled hexachlorobiphenyl (HCBP) to benthic amphipods from 
experimentally contaminated natural sediments. They found that organisms 
exposed to water and contaminated sediments had higher (from 2.3 to 10 
times) concentrations than did organisms exposed only to water. 
Furthermore, they found that the ratio of total bioaccumulation to 
Inoconcentration decreased with exposure time. This implied that the 
importance of food uptake decreased with exposure duration, and that the 
ultimate equilibrium state of the animal is independent of the uptake 
pathway. 

Schnoor (1982) attempted to determine the applicability of 
laboratory-derived BCFs for predicting actual field BCFs for a variety of 
hydrophobic organic chemicals in fish. He concluded that field BCFs, when 
normalized on a lipid basis, were within a factor of 1 to 4 of 
laboratory-derived BCFs. He suggested that laboratory-derived BCFs could 
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Ik? used in fate and effects modeling with order of magnitude accuracy, or 

better. 

INTERPRETIVE FRAMEWORK FOR BENTBIC BIOACCUMULATION 

In attempting to relate sediment exposure concentrations of organic 
chemicals to concentrations in benthic organisms and fish, two questions 
immediately arise. The first question concerns the validity of extending 
existing relationships between BCF and K from pelagic animals to benthic 
animals. The second question concerns the characterization of phase 
partitioning of organic chemicals in the sediment. This issue arises 
because actual field data for interstitial water concentrations are 
virtually non-existent. 

Following the arguments of Mackay (1982) that bioconcentration is a 
simple phase partitioning phenomenon, it seems reasonable to apply 
existing relationships between BCF and K to benthic organisms and fish. 
The principal difficulty occurs in specification of an appropriate 
exposure concentration of the chemical in the interstitial water. In 
principle, equilibrium partition relationships could be used to estimate 
interstitial water concentrations from total bulk sediment concentrations 
and sediment solids concentration. However, the validity of conventional 
partitioning theory has been challenged and the entire area of phase 
partitioning of organic chemicals in bottom sediments is an unresolved 
research issue. This is discussed below in more detail. As a first 
approximation, the conventional theory of reversible equilibrium 
partitioning is assumed to be valid in the sediments. 

Connor (1984) has suggested that chemical concentrations in benthic 
fish should be referenced to sediment concentrations instead of directly 
to aqueous phase concentrations. This approach assumes that fish and 
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sediments will be in equilibrium through their individual equilibria with 
the aqueous phase. This assumption is consistent with the theoretical 
arguments by Mackay (1982) on multi-phase equilibrium partitioning. 

Connor originally developed a relationship for the fish/sediment ratio 
of a chemical in terms of K and the fraction organic carbon (f QC ) in the 
sediment. This relationship was based on a combination of two regression 
equations developed by Kenaga and Goring (1980) for BCF and sediment 
partition coefficient (K ,). Connor (1985) proposed that regression 
equations developed by Mackay (1982) for fish BCF, and by Karickhoff et 
al. (1979) for organic chemical phase partitioning, be combined to yield 
the following expression for fish/sediment ratio: 

BCF 0.076 m 

V = — E UJ 

h sed r oc 

Oimensionally, the ratio BCF/K , is equivalent to the ratio of chemical 
concentration in fish (wet weight) to chemical concentration in sediment 
(dry weight). This equation has appeal because of its simplicity and 
because it contains no explicit dependence on K or on interstitial water 

KJYI 

concentration. 

An important assumption in the Mackay-Karickhoff equation is that 
phase partitioning of organic chemicals in the sediment follows the 
relationship developed by Karickhoff et al. (1979) for more dilute 
suspensions. In particular, potential complications due to the effects of 
solids concentrations on equilibrium partition coefficients and to 
colloidal phases are not explicitly addressed. These are discussed below 
in more detail. 

The predictions of animal/sediment concentration ratio from the 
Mackay-Karickhoff equation are illustrated in Fig. 5.4. Consistent with 
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the definition of bioconcentration, and with the assumptions of the two 
component regression equations, this relationship provides a conservative 
estimate of actual environmental residues in benthic organisms and fish. 
Furthermore, it provides an approximate criterion for determining whether 
biomagnif ication may be contributing to observed accumulation levels. 
FIELD DATA FOR BENTHIC BIO ACCUMULATION 

This section contains an examination of available field data for 
animal/sediment concentration ratios of organic chemicals in benthic 
organisms and fish. The objectives of this analysis are to determine if 
actual residue values are consistent with the simplified theory underlying 
the Mackay-Karickhoff equation and to determine whether biomagnif ication 
might be occurring. 

Field data for chlorobenzenes, PCBs and hexachlorobutadiene (HCBD) in 
Lake Ontario oligochaetes and amphipods have been reported by Fox et al. 
(1983) (Figs. 5.5 and 5.6). The oligochaetes were primarily Tubifex 
tubifex and the amphipods were primarily Pontoporeia hoyi. Reported 
animal concentration data were converted from dry to wet weight by 
assuming that dry weight was 157. of wet weight (Oliver 1984) . The curve 
on each plot represents the predictions of the Mackay-Karickhoff equation. 

In general, the higher ratios correspond to the more highly 
chlorinated chemicals. In particular, the highest values in each case 
correspond to hexachlorobenzene (HCB) and PCB. These have the highest K 
values of the reported chemicals and they do not degrade significantly in 
the environment. The lower chlorobenzenes are thought to be metabolized 
by the benthic organisms and/or microbially degraded. 

Although the amphipods tended to have higher ratios than the 
oligochaetes, results for both groups of organisms were reasonably 
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consistent with predictions from the simplified partitioning theory. 
Mackay (1982) reported a standard error of 1.8 BCF in his regression of 

BCF vs K . Karickhoff et al. (1979) did not report the standard error 

o 
for their regression of sediment partition coefficient vs K Qw (r = 0.96). 

None of the oligochaete ratios, and only one of the amphipod ratios, 

exceed the Mackay-Karickhoff predictions by greater than a factor of 2. 

It is possible that differences in feeding habits are responsible for 
the observed differences in chemical accumulation. Oligochaetes live 
within the surficial sediment and continuously ingest sediment particles. 
Amphipods are found at the sediment-water interface and deposit feed 
directly on detritus. 

Data for PAHs in Lake Erie oligochaete worms and chironomid midges 
have been reported by Eadie et al. (1983b) (Figs. 5.7 and 5.8). The worms 
were primarily Limnodrilus and the midges were primarily Chironomus and 
Procladius. In general, the ratios for the worms and midges are 
consistent with the simplified theory. None of the worm ratios, and only 
two of the midge ratios, exceed the Mackay-Karickhoff predictions by 
greater than a factor of 2. 

Data for PAHs in Lake Michigan Pontoporeia have also been reported by 

Eadie et al. (1983b) (Fig. 5.9). These ratios appear somewhat higher than 

the predictions of the simplified theory. Two of the four data points for 

f = 0.004, and all six of the data points for f = 0.017, exceed the 

Mackay-Karickhoff predictions by greater than a factor of 2. It is 

< 
possible that biomagnif ication, as well as bioconcentration, are both 

contributing to the observed accumulation levels. 

Data for PCBs in benthic fish have been summarized by Connor (1984) 

and are re-plotted in Fig. 5.10. These data are from ten separate studies 
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and they include a laboratory aquarium, six marine sites and three 
freshwater sites (one in Lake St. Clair and one in Lake Erie). The 
organic carbon content of the associated sediments was not reported in any 
of the studies. The curve on the plot represents the predictions of the 
Mackay-Karickhoff equation for f = 0.01 (1 percent). These predictions 
are a reasonable upper bound from the simple theory because the actual f 
values of the sediments were probably greater than 0.01, and because 
liigher values of f would correspond to lower predictions for the 
fish/sediment ratio. 

Connor plotted the data as a function of flushing time because he 
postulated that a certain minimum water residence time might be required 
to ensure equilibrium conditions between sediments and fish. He 
speculated that the fish/sediment ratios appeared to plateau at flushing 
times greater than 100 days, suggesting that the time necessary for 
fish-sediment equilibrium in the field may be of the same order as in 
laboratory bioconcentration studies. 

Apart from any possible relationship between fish-sediment 
concentration ratios and flushing time, there appears to be much scatter 
in the data. Furthermore, four of the ratios exceed the upper bound 
prediction of the Mackay-Karickhoff equation by greater than a factor of 
2, and one of them exceeds this prediction by an order of magnitude. If 
this simplified partitioning theory is valid then it appears that 
biomagnification may account for a significant portion of the total 
accumulation of PCBs by benthic fish. Clearly, the limited data presented 
here are insufficient to resolve this important issue. 
RELATIONSHIP BETVEEN BI0ACCUMULATI0N AND SEDIMENT EXPOSURES 

Two potential weaknesses in the theory underlying the simplified 
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Mackay-Karickhoff equation are lack of explicit dependence of K , on 
solids concentration and lack of consideration of potential colloidal 
phases. O'Connor and Connolly (1980) synthesized results from a large 
number of independent laboratory studies and demonstrated that an inverse 
relationship existed between equilibrium partition coefficients and the 
concentrations of adsorbing solids. Horzempa and DiToro (1983) found that 
although variations in solution chemical composition and kinetic effects 
influenced partitioning, neither factor appeared to completely account for 
the magnitude of the observed sediment concentration effect in their 
experiments. 

The potential influence of suspended solids concentrations on 
partition coefficients is not limited to laboratory experiments. Capel 
and Eisenreich (1985) found an inverse log relationship between PCB 
partition coefficients and suspended solids concentrations in Lake 
Superior. This behavior was partially attributed to differences in 
organic carbon content of the particles in different regions of the lake 
and to artifacts in the filtration process used. Swackhamer and Armstrong 
(1987) found a similar relationship in Lake Michigan. 

Some workers have explained the influence of solids concentration on 
partition coefficients in terms of "third phases". Voice et al. (1983) 
studied the sorption of four hydrophobic chemicals on three different 
sediments from Lake Michigan. They found significant increases in 
equilibrium partition coefficents as solids concentration decreased. They 
suggested that this phenomenon was caused by the presence of 
microparticulates contributed by the adsorbing solids and not removed from 
suspension in the separation procedure. Gschwend and Wu (1985) studied 
the sorption of several PCB isomers on bottom sediments from Lake Superior 
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and the Missouri River. They contended that if precautions are taken to 
eliminate or account for non-settling (or non-filterable) 
microparticulates or organic macromolecules which remain in the aqueous 
phase during laboratory experiments, then equilibrium partition 
coefficients will remain constant over a wide range in suspended solids 
concentration. 

The potential association of hydrophobic organic chemicals with 
dissolved organic matter can also confound the interpretation of 
experimental sorption data. Hassett and Anderson (1979) reported that 
cholesterol could be bound by dissolved organic matter (DOM) and 
stabilized in the aqueous phase. Chiou et al. (1986) found significant 
enchancements in solubility due to association with DOM of soil and 
aquatic origins. They contended that this phenomenon could cause the 
apparent dependence of partition coefficients on sediment concentrations. 
Carter and Suffet (1982) suggested that a significant fraction of DDT in 
natural waters may be bound to dissolved humic material, and that this 
could cause an apparent particle concentration effect on equilibrium 
partition coefficients. Landrum et al. (1984) showed that the equilibrium 
partition coefficient for binding of a variety of hydrophobic chemicals to 
humic acid and dissolved organic carbon (DOC) in natural waters was 
inversely proportional to the concentration of DOC in solution. 

The presence of colloidal phases may significantly affect the 
bioavailability of organic chemicals. Landrum et al. (1985) suggested 
that reduction in bioavailability in the presence of DOC may explain 
anomalies observed in predicting bioaccumulation of chemicals with log K 
> 6 (e.g. Oliver 1984). He postulated that much higher proportions of 
high log K chemicals would be complexed by DOC as compared to lower K 
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chemicals. Complexation by DOC would tend to reduce the true aqueous 
phase concentration of the chemical. Landrum et al. (1987) demonstrated 
that DOM in actual interstitial waters from Lake Michigan sediments 
reduced the bioavailability of PAH and PCB to Pontoporeia. 

If partition coefficients are, in fact, lower in the sediment than 
predicted by conventional partitioning theory (e.g., Karickhoff et al. 
1979) , then benthic organisms would be exposed to higher aqueous phase 
concentrations in the interstitial water than in the water column. 
Conversely, if partitioning to DOM/DOC is significant in the sediment, 
then higher interstitial water concentrations may not represent true 
aqueous phase concentrations. If only the true aqueous phase 
concentrations of organic chemicals are bioavailable, then these two 
phenomena may tend to oppose each other. 

DiToro (1985) has proposed a mechanistic particle interaction model 
for the reversible component partition coefficient. He successfully 
applied this model to data for a large number of neutral hydrophobic 
organic chemicals. Figure 5.11 contains a comparison between the 
predictions of this model for PCB partition coefficient, and data for the 
water column (Eadie et al. 1983a) and sediments (Rice et al. 1983) in Lake 
Michigan. The data indicate that apparent partition coefficients are 
significantly lower in the sediments than in the water column. The model 
predictions agree with the data for the water column, but they are 
substantially lower than the data for the sediments. It is not possible 
to determine if the measured sediment partition coefficients were 
underestimated due to the presence of a non-separable third phase. 

It is instructive to examine the relationship between concentrations 
of organic chemicals in benthic organisms and the corresponding 



equilibrium aqueous phase concentrations. This approach has been used in 
the water column to estimate aqueous phase concentrations of DDT and PCBs 
from fish body burdens (Bierman and Swain 1982; Rodgers and Swain 1983). 
Figure 5.12 illustrates the hypothetical animal-water relationships for 
PCBs (log BCF =5.15 from Mackay 1982) using the oligochaete and amphipod 
data from Lake Ontario, and the interstitial water concentration data from 
Lake Michigan. 

Assuming that bioconcentration is the only process responsible for 
accumulation in the animals, results indicate that the oligochaetes and 
amphipods would be in equilibrium with aqueous phase concentrations of 1 
to 20 ng/L. Conversely, under the same assumption, the interstitial pore 
water concentrations would result in animal concentrations of 10 to 100 
ppm (wet weight) at equilibrium. Water column concentrations of PCBs in 
the Great Lakes range between 1 and 10 ng/L (e.g., Rice et al. 1983). 
Animal concentrations of PCBs generally do not exceed 10 ppm (wet weight) 
(Great Lakes Water Quality Board 1983). 

These results imply that true aqueous concentrations in the 
interstitial water are not substantially different than aqueous 
concentrations in the water column, or that benthic organisms are in 
equilibrium with water column concentrations instead of sediment 
concentrations. Alternatively, if benthic organisms are at equilibrium 
with reported interstitial water concentrations, then these concentrations 
are not completely bioavailable. 
UNRESOLVED RESEARCH ISSUES 

There are many unresolved research issues in the area of organic 
chemical accumulation by benthic organisms and fish. Historically, much 
more emphasis has been placed on the pelagic food chain. Some of these 
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unresolved issues relate to the sediment exposure environment and others 
relate to the linkage between sediment exposures and bioaccumulation. 
Still other issues relate to the very structure of benthic food chains. 
The major unresolved issues related to benthic food chains are: 

1. The number and types of trophic linkages within the benthic food 
chain . 

2. The number and types of linkages between the benthic and pelagic 
food chains. 

3. Identification of the major carbon transport pathways through the 
benthic food chain. 

4. Identification of the minimum number of feeding types and trophic 
compartments required to account for the majority of the carbon 
transport through the benthic food chain. 

The major unresolved issues related to the sediment exposure 
environment are: 

1. Characterization of chemical phase partitioning. Methods are 
needed for estimating interstitial water concentrations from bulk 
sediment properties that can be more easily measured. 

2. Determination of the significance of colloidal phases and 
complexation to macromolecules. 

3. Characterization of the bioavailable component of total sediment 
exposure concentrations. 

4. Quantification of the physical-chemical alterations to the 
sediments by benthic organisms. 

5. Determination of the significance of microbially-mediated 
degradation of organic chemicals in sediments. 

Some of the major unresolved issues related to bioaccumulation by 
benthic organisms and benthic fish are: 

1. Validity and limits of applicability of existing empirical 
correlations between bioconcentration factors and chemical 
exposure concentrations. 

2. Potential improvements to simple empirical correlations by 
normalizing to lipid content or to other characteristic biological 
parameters. 
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3. Determination of ingestion, assimilation, metabolism and excretion 
rates of organic chemicals by benthic organisms and fish. 

The single principal obstacle to understanding and describing organic 
chemical dynamics in benthic food chains is lack of adequate field data 
for sediment exposures and animal accumulation levels. Field data are 
required in order to make initial judgements on the appropriate degree of 
model complexity for benthic bioaccumulation and potential coupling to the 
pelagic food chain. Ultimately, field data provide the final test for 
different hypotheses on model structure and process mechanisms. 
SUMMARY AND CONCLUSIONS 

There is considerable uncertainty regarding the structure and function 
of the benthic food chains in the Great Lakes. Furthermore, potential 
coupling between the benthic and pelagic food chains is poorly understood. 

Benthic organisms can have significant effects on the physical and 
chemical properties of sediments and on constituent mass fluxes across the 
sediment-water interface. Models for microscopic scale processes of 
sorption kinetics and equilibrium, and for exposure of bed particles to 
overlying waters are only in the early stages of development. Sediment 
pelletization may mitigate sediment mass fluxes by trapping sorbed 
chemicals and significantly retarding chemical release into the aqueous 
phase at the sediment-water interface. 

Environmental contaminants can be accumulated by benthic organisms and 
fish through direct uptake from water and through food chain uptake. The 
primary mode of accumulation in a particular case depends on the 
properties of the chemical, the relative concentrations of the chemical in 
the water and food, and on the trophic position of the animal. Under 
certain conditions, octanol-water partition coefficients can be good 
predictors of bioconcentration factors in fish and benthic organisms. 
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These factors are only conservative estimates of actual environmental 
accumulation levels because these levels are the result of potential food 
uptake, as well as bioconcentration. 

A simplified partitioning relationship was used to examine the 
available field data for animal/sediment concentration ratios of organic 
chemicals in benthic organisms and fish. Reported ratios for oligochaete 
worms in Lakes Ontario and Erie, and chironomid midges in Lake Erie, were 
consistent with biococentration due to simple phase partitioning. 
Reported ratios for amphipods in Lake Ontario were also consistent with 
simple phase partitioning, however, ratios for amphipods in Lake Michigan 
appeared to indicate that biomagnification may be important. Data for 
benthic fish were extremely limited, however, they appeared to indicate 
that biomagnification accounted for a significant portion of total PCB 
body burdens. 

These preliminary results indicate that accumulation of organic 
chemicals in macrobenthic invertebrates may be approximated using 
empirical relationships for bioconcentration. However, the description of 
accumulation levels in benthic fish may require a mechanistic food chain 
modeling approach. 

The possible dependence of phase partition coefficients on solids 
concentration and potential interferences from colloidal phases are major 
obstacles to the characterization of sediment exposures. The available 
data are extremely limited, however, they appear to indicate that 
partition coefficients for PCBs are significantly lower in the sediment 
than in the water column. This implies that benthic organisms are exposed 
to much higher aqueous phase concentrations in the interstitial water than 
pelagic organisms exposed to the water column. However, reported values 
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for organic chemical accumulation in benthic organisms are not consistent 
with substantially higher exposures, thus indicating that interstitial 
water concentrations may not be completely bioavailable. 

The single major obstacle to understanding and describing organic 
chemical dynamics in benthic food chains is lack of adequate field data 
for sediment exposures and animal accumulation levels. 
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FIGURE CAPTIONS 

Fig. 5.1. Hypothesized trophic links in the Lake Ontario deepwater food 
web. Percentages suggest proportional contribution of food 
from various sources to trophic compartment (Flint, 1986). 

Fig. 5.2. Theoretical model of an annual carbon production and energy 
flow food web for the deepwater habitat of Lake Ontario. All 
compartment production and carbon flow values are represented 

in g C/m 2 /yr (Flint, 1986). 

Fig. 5.3. Simplified food chain for Lake Ontario (Oliver and Niimi, 
1987). 

Fig. 5.4. Predicted relationship between animal/sediment concentration 
ratio for organic chemicals and fraction sediment organic 
carbon using the Mackay-Karickhoff equation (see text). 

Fig. 5.5. Organochlorine contaminants in Lake Ontario oligochaetes. 

data from Fox et al. (1983). Curve is predicted 

animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 

Fig. 5.6. Organochlorine contaminants in Lake Ontario amphipods. Data 
from Fox et al. f 1983) . Curve is predicted animal/sediment 
ratio from Mackay-Karickhoff equation (see text). 

Fig. 5.7. Polyaromatic hydrocarbons (PAHs) in Lake Erie worms. Data 

from Eadie et al. (1983b). Curve is predicted 

animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 

Fig. 5.8. Polyaromatic hydrocarbons (PAHs) in Lake Erie midges. Data 
from Eadie et al. (1983b). Curve is predicted 
animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 

Fig. 5.9. Polyaromatic hydrocarbons (PAHs) in Lake Michigan 

Pontopereia. Data from Eadie et al. (1983b). Curve is 

predicted animal/sediment ratio from Mackay-Karickhoff 
equation (see text) . 

Fig. 5.10. Total polychlorinated biphenyls (PCBs) in benthic fish. Data 
from Connor (1984). Curve is predicted animal/sediment ratio 
from Mackay-Karickhoff equation (see text) for 17. sediment 
organic carbon. 

Fig. 5.11. Partition coefficient for PCBs as a function of solids 
Concentration. Data for total PCBs in water column (W) for 
Lake Michigan from Rice et al. (1983); Data for total PCBs in 
sediment (S) for Lake Michigan from Eadie et al (1983a). 
Curve is predicted partion coefficient for hexachlorobiphenyl 
from DiToro (1985). 
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Fig. 5.12. Hypothetical animal-water equilibrium relationships for total 
PCBs. Log BCF = 5.15 from Hackay (1982). Oligochaete (0) 
and amphipoa (A) PCB concentrations from Fox et al. (1983) 
for Lake Ontario. Pore Water PCB concentrations (P) from 
Eadie et al. (1983a) for Lake Michigan. 
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FIGURE 5.1: Hypothesized trophic links in the Lake Ontario 
deepwater food web. Percentages suggest proportional 
contribution of food from various sources to trophic com- 
ponent. (Flint, 1986) 
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FIGURE 5.2: Theoretical model of an annual carbon production 
and 'energy flow food web for the deepwater habitat of Lake 
Ontario. All compartment production and carbon flow values 
are represented in g C/m2/yr. 
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Simplified food chain for Lake Ontario (Oliver and Hilml 1987) 
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Fig. 5.4. Predicted relationship between animal/sediment concentration 
ratio for organic chemicals and fraction sediment organic 
carbon using the Mackay-Karickhoff equation (see text). 
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Fig. 5.5. Organochlorine contaminants in Lake Ontario oligochaetes. 
data from Fox et al. (1983). Curve is predicted 
animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 
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Fig. 5.6. 



Qrganochlorine contaminants in Lake Ontario amphipods. Data 
from Fox et al. (1983). Curve is predicted animal/sediment 
ratio from Mackay-KariCKhoff equation (see text). 
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PAHa in Laka Erie Worm* 
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Fig. 5.7. Polyaromatic hydrocarbons (PAHsJ in Lake Erie worms. Data 

from Eadie et al. (1983d). Curve is predicted 

animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 
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PAHs in Lake Eria Midgas 
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Fig. 5.8. 



Polyaromatic hydrocarbons (PAHs) in Lake Erie midges. Data 
from Eadie et al. (1983b). Curve is - predicted 
animal/sediment ratio from Mackay-Karickhoff equation (see 
text) . 
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PAHs in Lake Michigan Pontoparsia 
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Fig. 5.9. Polyaromatic hydrocarbons (PAHs) in Lake Michigan 

Poniopereia. Data from Eadie et al. (1983b). Curve is 

predicted animal/sediment ratio from Mackay-Karickhoff 
equation (see text). 
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Total PCBa in Benthic Fish 
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Fig. 5.10. Total polychlorinated biphenyls (PCBs) in benthic fish. Data 
from Connor (1984). Curve is predicted animal/sediment ratio 
from Mackay-Karickhoff equation (see text) for 17. sediment 
organic carbon. 
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PCB Partition Coefficient as a 
Function of Solids Concentration 
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Fig. 5.11 



Partition coefficient for PCBs as a function of solids 
Concentration. Data for total PCBs in water column (V) for 
Lake Michigan from Rice et al. (1983); Data for total PCBs in 
sediment (S) for Lake Michigan from Eadie et al (1983a). 
Curve is predicted pr.rtion coefficient for hexachlorobiphenyl 
from DiToro (1985). 
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Hupat.heti.cal Animal -Water Equilibrium 
Relationship* far Total PCBs 
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12. Hypothetical 
PCBs. Los 



aniniai-water equilibrium relationships for total 
BCF = 5.15 from Mackay (19S2). Oligochaete (0) 
and amphipod (A) PCB concentrations from Fox et al. (1983) 
for Lake Ontario. Pore Water PCB concentrations (P) from 
Eadie et al. (19S3a) for Lake Michigan. 
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CHAPTER 6 



FEEDING RELATIONSHIPS AND CONTAMINANT PATHWAYS PART I 



G.D. HAFFNER 

GREAT LAKES INSTITUTE 
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INTRODUCTION 

To consider feeding relationships in biological systems one must 
recognize that biological systems are often built upon the concepts that 
energy must be transferred efficiently from one trophic level to another, 
and that the systems must also be stable over time and space. Long food 
chains are inherently unstable. They represent the most efficient pathway 
to transfer energy and therefore contaminants from the primary producers 
to the top predator. Simple food chains however require such close 
synchrony that any unusual event can interrupt the energy flow such that 
the resources at the upper level of the food chain can be seriously 
depleted. Resource managers in the past have given us prime examples as 
to what happens to ecological systems when food chain models have been 
selected to manage natural resources. Fisheries at one time spoke of 
maximum sustainable yield, and of course now speak of restoration and 
rehabilitation. The agricultural business selected plants for their fast 
growing abilities not for their resistance to disease or insects, thus our 
ever growing dependence on pesticides and fungicides. To consider feeding 
relationships one must consider food web models, and the positions of 
species within the food web. These are naturally much more complex, but 
still a great deal can be gained by looking at food webs as an integral 
part of environmental resource management. 
FOOD VEBS 

Let us first of all look at a generalized food web, and the properties 
one must consider of a food web. The first property is that the number of 
trophic levels that can be sustained in a food web are not a simple 
function of total energy flow. Indeed there appears to be a paradox in 
that a diversity of species is essential to develop a stable, multi-tiered 
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food web. This diversity of species obviously has an energy cost to 

maintain the interactive and interdependent relationships of species that 

comprise the food web. Energy spent in this manner is not available at 

other trophic levels, thus there is an energy cost for stability in food 

webs. 

Stability in ecological systems is probably one of the most 

controversial areas one can enter, yet one cannot discuss ecological 

systems without recognizing there is a stability, i.e. we do recognize 

pattern in ecological systems. Some researchers such as Cairns and Van 

der Schalie consider there to be an "ecosystem assimilation" potential for 

dealing with (i.e. tolerating) the direct effects of toxic substances. To 

avoid some of the confusion associated with the word "stability" it is 

often better to consider two separate attributes of stability. The first 

is "perseverance" or the fact that food webs and ecosystems tend to 

conserve the number of species of which they are composed. The maximum 

number of potential links amongst all the species is given by 
n (n - 1) 

One of the better reviews of food webs with respect to perseverance is 
that of Briand (1980) who studied numerous terrestrial and aquatic food 
webs and concluded one could possibly develop a matrix to illustrate the 
connectance of species within a food web. It is probably not wise to 
generalize that the species with the greater number of connections are 
more important than those with less, but the loss of species with the 
greater number if interactions can possibly have a more profound effect on 
the long term stability of the food web. Should such species be 
physiologically more susceptible to the exposure of toxic substances, one 
would have to question what Cairns and Van der Schalie meant when 
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discussing "ecosystem assimilation capabilities". The loss of species to 
a food web undermines the stability of that food web, and there are good 
ecological reasons why environmental managers should return to the old 
concept of protecting the most sensitive species. 

The second property of stability is probably well recognized by those 
who frequently discuss "food chain" dynamics, for it is closely related to 
the transfer of energy from one trophic level to another. This is the 
"constancy" of populations within an ecodeme. It is not surprising to see 
the number of discussions being held today on the Great Lakes ecosystem as 
to whether food chains are controlled from the top down or from the bottom 
up. I would like to propose today a comprise between food webbers and 
food chainers but unfortunately, ecology does not lend itself to 
compromise. Constancy and perseverance tend to result in very different 
management options. It is quite different to regulate or maximize the 
abundance of key economic species, than to try to maintain the diversity 
of the food web to which that species belongs. 

Another important aspect not captured in the food chain concept is 
that there are many different associations that can exist among species. 
For example in Fig. 6.1 we can see that predation is only part of the 
species interactions found within a food web, and that these relationshps 
can have profound effects on energy flow or material transfer within an 
ecodeme. Other relationships such as competition, amensalism, parasitism, 
neutralism, commensal ism, and mutualism should also be considered as such 
interactions help determine the structure and function of food webs, yet 
are ignored within the discussion of simple food chains. 
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FOOD VEB MODELLING 

When attempting to model ecosystems, one is essentially trying to 
determine energy flows, recognizing that although we often draw energy 
"loops", energy is continually lost as heat, friction and advective 
transport of material out of the food web (e.g. loss of organic carbon to 
the deep sediments) . Therefore when one is confronted with food web 
diagrams or models, one must first of all realize that what is being 
modelled is often not species interactions or interdependencies, but is 
essentially an energy budget. In order to model a representative energy 
budget for a food web one has to have measurements of the rate limiting 
steps illustrated in Fig. 6.2. 

Factors that control the size and shape of the energy budget are: 

P 
Ecological efficiency = p — — j- 

x 

where P is the production at trophic level x 

I x ' 

Exploitation efficiency = p — — j- 

where I is the ingested portion at trophic level x 

A 

Assimilation efficiency = t — — j- 

where A is the assimilated portion of food ingested 

A 

P x 

Net production efficiency = t- 

x 

Gross production = I x - E x 

where E is the excretory loss at trophic level x 

Net production = G.P. (gross production) minus R x (respiratory loss) at 

trophic level x. 

One cannot look at such a list and not observe how important the word 
"efficiency" is. For example, assimilation efficiency varies in the 
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literature for different species from 157. to 90%. As there is a finite 
amount of energy entering the system and ecological efficiency is low, the 
number of trophic levels in a food web are limited. As mentioned earlier, 
ecological efficiency must compensate somewhat for system stability. 

This limited number of trophic levels gave rise to the concept of the 
Eltonian pyramid, yet these can be quite different in terrestrial and 
aquatic systems Fig. 6.3. The production biomass ratio for terrestrial 
systems is much smaller than that for aquatic systems, and the length of 
the subsequent "food chain" in aquatic systems is frequently longer than 
that in terrestrial ecosystems. 

The problem with Fig. 6.3, however, is that it tends to imply that 
energy can move only to the adjacent trophic level (steps cannot be 
skipped). Often however, as steps are skipped, and combined with the 
transfer of energy by detritivors into various tiers, one can see why 
ecological efficiency is very difficult to measure. This is an obvious 
area requiring more research if models of energy or contaminants can meet 
Levins (1960) criteria of being precise, general and realistic. Whether a 
heuristic model with driving processes as shown in Fig. 6.2 can simulate 
ecological systems in a realistic and precise manner remains to be seen. 
Interestingly, despite the need for such knowledge, this tends to be one 
of the poorest funded research areas simply because one must produce a 
monograph as opposed to a series of related papers. 
FOOD WEB DYNAMICS AND CONTAMINANTS 

The factors regulating food web models noted in Fig. 6.2, are 
essential for modelling contaminants, partitioning and concentrations in 
the food web. Contaminant burdens are measured often with regard to the 
ratio of contaminant mass to biological mass. If there is rapid recycling 
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of the biological mass in aquatic systems as observed in Fig. 6.3, and as 
the rate of recycling is decreased from one trophic level to the next 
there is an obvious problem in representing contaminant information in 
such a static fashion. An exposure at higher trophic levels can result in 
higher body burdens than at lower trophic levels. Simple contaminant 
bioaccumulation models often ignore this fact, and may result in false 
estimates of contaminant transfer. 

Other important aspects of contaminant transfer is prey selection, 
seasonality and life cycle changes. It is not surprising to observe that 
aquatic contaminant transfers are often regulated at this level of the 
phytoplankton. Despite very Tow body burdens, any mass balance model 
based on production and not standing stock would greatly effect the 
prediction of many food chain models (Flint 1986). 
CONCLUSION 

I would hesitate to recommend we try to model and monitor all 
compartments and the species within food webs (indeed, during their life 
cycle, many species can change compartments). I think we must start with 
a simple model before we can start making testable hypotheses as to the 
movement and accumulation of contaminants in food webs. Perhaps food webs 
or food chains are the wrong concept to be using especially for the Great 
Lakes ecosystem until some basic and long overdue research on primary 
productivity, zooplankton grazing and fish predation has been done. We 
have ignored these factors during the eutrophication issue only to end up 
having "top down or bottom up" food web workshops. It would appear that, 
unless environmental managers wish to end up in the same position again, 
we must finally start doing some basic ecological research in the Great 
Lakes ecosystem. 
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FIGURE CAPTIONS 

Fig. 6.1 



Fig. 6.2. 



Fig. 6.3. 



A summary of biological interactions that can regulate the 
relative abundance and composition of communities. 

Factors regulating energy flow, carbon flow and potentially 
contaminant flow in ecosystems. 

Production in trophic level x 

Amount ingested 

Amount of energy assimilated 

Amount of energy excreted 

Amount of energy lost through respiration 

Next lower trophic level 



x + 1: 



Energy flow models of terrestrial and aquatic ecosystems. 
Note the different production: biomass ratios as discussed 
in the text. 
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BIOLOGICAL RELATIONSHIPS 



species no. 1 



mutualism 



CM 

8 

o 

CD 

or 



o commensalism neutralism 



predation 
- parasitism ammensalism competition 

Figure 6.11 A summary of biological interactions that can regulate the 
relative abundance and composition of communities. 



129 



REGULATING FACTORS 
Ecological Efficiency = Px 

Px+I 

Exploitation Efficiency = lx 

Px~T 

Assimilation Efficiency = Ax 

lx 

Net Production Efficiency a Px 

Ax 

Gross Production = / x - Ex 

Afef Production = Ax-Rx 

= G.P-Rx 



Fig. 6.2. Factors regulating energy flow, carbon flow and potentially 
contaminant flow in ecosystems. 



P 
x 



Production in trophic level x 
Amount ingested 
Amount of energy assimilated 
Amount of energy excreted 
Amount of energy lost through respiration 
x + 1: Next lower trophic level 
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ENERGY FLOW MODELS 



TERRESTRIAL 



OWLS 
FOXES 
RABBITS 
PLANTS 



Prod/Biomass - 0.33 



AQUATIC 



m- . * 



PISCIVOROUS FISH 



PLANKTIVOROUS FISH 



CARNIVOROUS ZOOPLANKTON 



HERBIVOROUS ZOOPLANKTONfl 



ALGAE 



Prod/ Bio mass -25-42 



Fig. 6.3. Energy flow models of terrestrial and aquatic ecosystems. 

Note the different production: biomass ratios as discussed 
in the text. 



131 



CHAPTER 7 

FEEDING RELATIONSHIPS PART II: 
THE USE OF ECOLOGICAL FOOD WEB MODELS 



I. MICHAEL WEIS 

GREAT LAKES INSTITUTE 

AND 

DEPARTMENT OF BIOLOGICAL SCIENCES 

UNIVERSITY OF WINDSOR 

WINDSOR, ONTARIO, CANADA 

N9B 3P4 



INTRODUCTION 

Ecologists can make many contributions to the development of 
contaminant food web models from points of view other than those of 
mathematical modellers or fisheries biologists. One of those roles is as 
gadflies, ensuring that components of the ecosystem are all considered in 
the construction of models, and that the effects of simplifications 
necessary to construct mathematically tractable models are also 
considered. The object of model development is to permit predictions of 
the effects of changes in the ecosystem or the forcing functions on 
particular components of the system. The objective of food web models in 
this case, at least in part, is to permit prediction of contaminant 
burdens in sport fish such as pickerel, pike and bass. There are many 
models which could be employed to make those predictions, ranging from 
estimates based only on contaminant levels in sediments and water, through 
various kinds of simple, linear food chain models, to arbitrarily complex 
models representing the best available models of interlinked food webs. 

Simple models have considerable appeal to those who must construct and 
test them. More complex models appeal to biologists, who typically 
believe they better represent interactions in the real world. The usual 
systems, or box-and-arrow model of a food web represents a compromise 
between biologically complex reality and computational tractability. To 
ensure that those who construct box-and-arrow models recognize that these 
models are usually simplifications of biological reality, a few examples 
of evidence of the complexity in the biology of one of the target taxa 
will serve. A survey of the published literature over the 1984 to 1986 
period studying the ecology of large and small mouth bass ( Micropterus 
salmoides and dolomieui respectively) indicates a number of ecological and 
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behavioural complexities difficult to include in a box-and-arrow model of 
the sort proposed. 
SYSTEM COMPLEXITIES 

One kind of complexity results from temporal change in the ecology of 
bass. There is evidence of a consistent pattern of age and/or size 
dependent shift in the habitat selection of both large and small mouth 
bass. Each species moves from shallow to deeper water with increasing age 
(Vanjala et al. 1986, Rankin 1986), and chooses specific bottom types as 
well. Spatial shifts need not have an impact on the food web model, but 
in parallel with habitat shifts various diet shifts also occur. Some 
studies of bass diets include suggestions of a simple correlation between 
prey size and predator size, indicating shifts in relative abundances of 
prey species in the diet with predator growth, but continuing use of the 
same spectrum of prey species. For example, large mouth bass are reported 
to consume at least some gizzard shad even as young-of-the year, but. 
gizzard shad became the dominant component of the gut contents in later 
years in Lake Shelbyville, Illinois (Storck 1986). That diet shift could 
be accommodated in a box-and-arrow model through age-dependent variation 
in transfer rates. However, in other cases the diet shifts reported are 
more radical, and represent transfers from completely different boxes at 
different ages. For example, large mouth bass have been reported to shift 
from a completely invertebrate diet when young to a completely piscivorous 
diet as mature adults (Miller and Starck 1984) . This complexity could be 
incorporated by considering different age classes as separate components 
of the ecosystem, but would require determination of a large number of 
inter-compartmental transfer rates. 

Other complexities are related to spatial variation in the ecosystem 
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structure and the behavioural responses of bass to those changes. Large 
mouth bass increase the breadth of their diets in response to increases in 
the density of aquatic vegetation in a laboratory study (Anderson 1984). 
This relationship between vegetation and diet was argued to result from 
the optimization of foraging tactics. Quantitative change in diet breadth 
is difficult to monitor in the field. Qualitative change, for example, 
change in preference from one species to another can be apparent in the 
field. Schramm and Zale (1985) indicate that with a cover of aquatic 
vegetation large mouth bass prefer bluegills, but in the absence of that 
cover their preferred food was telapia. Obviously habitat effects will be 
site specific. However, since adult bass move over a considerable area, 
and thus a range of habitat conditions, a quantitative estimate of diet 
could be problematic, and result in errors in the predictions of a 
contaminant web model. 

It seems likely that with sufficient biological data the kinds of 
temporal and spatial complexity readily evident in the literature relating 
to the diets of bass will also be found in the relationships among 
components at all levels of the food web. Further, the details of these 
complications are likely to show site specific differences. While each of 
the complexities considered, and this is clearly an incomplete list, make 
the development of a realistic, quantified box- and- arrow model difficult, 
they should also make evident the inadequacy of simpler models in their 
attempts to make accurate predictions. 
SYSTEMS MODELS 

Have box-and-arrow models been successful? If the answer is limited 
to attempts to estimate contaminant burdens, I know of no examples. There 
have, however, been numerous studies developing both mathematical models 
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and natural evidence for biomagnif ication. The field data are 
surprisingly equivocal concerning the widespread occurrence of 
biomagnif ication among species and trophic levels. There is general 
agreement that there is evidence of significant biomagnif ication of PCB's 
in higher trophic levels of both terrestrial and aquatic ecosystems. 
Sparrow hawks had more than 50 times the concentrations of PCB's in 
tissues of insectivorous birds in England (Presst et al., 1970). Thomann 
(1981) found contaminant levels of the top predators in aquatic ecosystems 
were similarly the result of biomagnif ication. The results of studies of 
lower trophic levels are much less clear. While biomagnif ication is tacit 
evidence of the importance of diet in contaminant uptake, lower trophic 
levels may accumulate a significant portion of their contaminant burdens 
through direct exposure, e.g. through uptake from water passing over gill 
structures in aquatic invertebrates. Many studies do not (and, if direct 
exposure is the principal route of uptake, should not) show evidence of 
biomagnif ication in lower trophic levels (Nisbit and Sarofim 1972, 
Ellgehausen et al. 1980, Hunter et al. 1980, and Clayton et al. 1977). 
Taken together these papers have studied a broad range of types of aquatic 
ecosystems from streams to marine environments. 

As a result of the uncertainty concerning the relative importance of 
dietary and direct environmental exposure to contaminants, mathematical 
models of the movement of contaminants among components of an ecosystem 
generally include terms for both sources. Shaw and Connell (1986) review 
a number of approaches. These approaches, as seems generally the case in 
contaminant studies, consider food chains, not food webs. They are, 
therefore, inadequate to model the movement of contaminants through the 
frequently complex pathways of either terrestrial or aquatic food webs. 
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This is the point at which ecologists may be able to aid in studies of 
transfers among compartments in real ecosystems. 

In the context of ecological, energetic food webs there are many 
examples, and many ways in which these can assist in the development of 
contaminant models. One of the major problems in developing contaminant 
models is the estimation of transfer rates among compartments. When 
energetic transfer rates are known, two approaches become possible: 1) 
Important links in the system can be identified from energetic transfer 
rates, and the contaminant web model simplified on ecologically justified 
grounds by excluding unimportant links from the final model. 2) Where 
contaminant transfer rates may be both difficult and expensive to measure, 
coarse estimates can be derived from the energetic web. These estimates 
assume that the lipid content of prey organisms carries an equilibrated 
concentration of the contaminant under study. If the lipid biomass 
transferred between boxes (trophic levels, or whatever units are being 
used in the model) can be estimated, then the contaminant transferred can 
also be estimated. These estimates are derived from the biomass transfer 
rates (0 ) , the percent lipid content of the prey organism (1 ), and an 
estimate of the lipofelicity of the contaminant under study, possibly K . 
The transfer rate for the contaminant, estimated from the energetic food 
web, is then, 
T = x 1 x K (1) 

X X OW v ' 

This approach to estimation is coarse and incomplete. Transfer rates 
among living components of the ecosystem are, of course, only part of the 
story. A complete model must include various losses which occur in 
addition to transfers to higher trophic levels. These include natural 
death rates, rates of metabolic degradation of contaminants, and excretion 
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rates. Contaminant movements due to excretion and organism death could be 
approached by parallel analogies to energetic models. Metabolic 
degradation has no energetic analogy, and is a complex problem for any 
method of food web modelling. This method does not necessarily 
distinguish potentially important species or individual differences in 
assimilation efficiency. However, it has the advantage of providing 
estimates necessary to the development of a screening model with relative 
ease. The necessary parameters have, in many cases, already been measured 
and published. 

Finally, energetic box-and-arrow models have been tested and verified 
against field measurements. Where simple models have been poor predictors 
of reality, the complications necessary to improve predictions have been 
determined and incorporated. Those complications provide useful clues to 
the difficulties likely to be encountered in developing a useful 
predictive model for the movement of contaminants through the food web. 
To see how such models are constructed and corrected, I will review 
briefly one model of a small lake ecosystem. The basis of this model is 
Lindeman's studies of Cedar Bog Lake (Linkeman 1941, 1942), in which much 
of the basic approach and terminology of trophic energy was developed. 
This trophic structure was quantified in a form useful to the development 
of contaminant web models by Williams (1971). 

Williams developed the final model in three stages. The first stage 
reiterated the three-level trophic structure consisting of compartments 
for plants including all photosynthetic organisms, herbivores and 
carnivores. For each compartment 1) input, 2) transfer to the level 
above, 3) respiratory loss, and 4) unspecified other losses (including 
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principally death) are balanced to bring the system to equilibrium. Each 
rate of transfer out from a box is linearly related to the standing crop 
of that compartment. The simple linear model which resulted did not 
accurately reflect the dynamics of measured standing crops. Therefore, 
Williams constructed a more detailed, 10-compartment model of the trophic 
structure of Cedar Bog Lake. 

To quantify the standing crops and transfer rates for the 10 
compartment model of Cedar Bog Lake, Williams had to make many 
assumptions, particularly about benthic components and the active layer of 
the sediment. The compartment sizes and flow rates when the differential 
equations describing the system had been solved for equilibrium 
conditions, could be used to assess the relative importance of different 
pathways (see Fig. 7.1). The largest compartment by a factor of 
approximately 50 is benthic plants. By virtue of that size, benthic 
plants are the chief contributors to the organic fraction of the ooze. 
Through indirect links the benthic plants are, therefore, important 
contributors to higher trophic levels. Major pathways from the ooze lead 
to upper trophic levels, e.g. piscivorous fish, by way of zooplankton and 
browsing herbivores. Both of these compartments received their largest 
inputs from the ooze. Other features evident in this web include 
differences in the major predators of nannoplankton (zooplankton 
predominate) and net plankton (plankton predators predominate) . This 
additional complexity still did not accurately predict the dynamics of 
standing crops in the compartments. 

To accurately portray those dynamics a third model, much like the 
second, but incorporating non-linearities in the responses of some 
compartments (e.g. self-inhibition or competition in the dynamics of plant 
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and swimming predator compartments) , differences in the temperature 
dependence (Q-.q) of responses of some compartments and the incorporation 
of seasonal cycles in the solar energy input and in the responses of some 
compartments produced a model which fit observed dynamics reasonably well. 
CONCLUDING COMMENT 

What have these complications demonstrated? A predictive model of the 
trophic web for Cedar Bog Lake had to incorporate a variety of realistic 
biological corrections, both in terms of the number and relationship of 
compartments, and in the internal dynamics of those compoartments, i.e. 
self-inhibition, to accurately model the dynamics of the ecosystem as a 
whole. The assumption that simple serial models can adequately predict 
the dynamics of contaminants in aquatic ecosystems is fraught with danger. 
Realistic models are usually, though unfortunately, realistically complex. 
REFERENCES 

Anderson, 0., (1984). Ecology, 65:851. 

Clayton, J.R., Pavlou, S.P., and Brietner, N.F., (1977). Environ. Sci. 
Technol., 11:676. 

Ellgehausen, G., Guth, G.A., and Esser, H.O., (1980). Ecotoxicol. 
Environ. Safety, 4:134. 

Hunter, R.G., Randolph, J.C., Carroll, J.H., (1980). Environ. Pollut. B, 
1:233. 

Lindeman, R.L., (1941). Amer. Midi. Natur., 26:636. 

Lindeman, R.L., (1942). Ecology, 23:399-418. 

Miller, S.J., and Strock, T.W., (1984). Trans. Am. Fish. Soc. , 113:571. 

Nisbit, I.C.T., and Sarofim, A.F., (1972). Environ. Health Perspect. , 
1:21. 

Presst, I., Jeffries, D.J., and Moore, N.W., (1970). Environ. Pollut., 
1:3. 

Rankin, E.T. , (1986). Trans. Am. Fish. Soc, 115:322-334. 



140 



Schramm, H.L., and Zale, A.V., (1985). Trans. Am. Fish. Soc, 
114:725-731. 

Shaw, G.R. , Connell. D.V., (1986). PCBs and the Environment, Vol. 1, 
J.S. Vaid (ed.), CRC Press, Boca Raton, Fl. } 135-142. 

Storck, T.V., (1986). Trans. Am. Fish. Soc, 115:21-27. 

Thomann, R.V., (1981). Can. J. Fish. Aquat. Sci., 38:280. 



141 



FIGURE CAPTIONS 

Fig. 7.1. Food web energy flows. 

Fig. 7.2. Major compartments and pathways in the Cedar Bog Lake 
ecosystem. For compartment sizes and non- linear equilibrium 
flow rates see Williams (1971). 
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FIG. 7-1 
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FIG* 7-2 Major compartments and pathways in the Cedar Bog Lake 
ecosystem. For compartment sizes and non-linear 
equilibrium flow rates see Williams (1971) . 
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SUMMARY OF FINDINGS 



The synopsis of the deliberation of the work groups as reported by the 

respective coordinators is attached in Appendix B. The major findings 

are: 

- The chemical classification based on K (the octanol water 

partition coefficient) can be used to assist in the selection of 
the appropriate model for estimating contaminant levels in biota. 

FOOD CHAIN/CHEMICAL CLASSIFICATION 



Log K 



ow 



Suggested Model 



Reference 



< 5 



5-7 



> 7 



Only phase partitioning required. 
Lipid normalization useful. 
Suggested screening model 
v - (K) (lipid fraction) c. 



R. Thomann (Chap. 3) 
D. Mackay (Chap. 4) 



Food chain model needed. 
Generic models are useful for screening. 
Site- and species-specific models are very 
data-intensive and are in the initial 
stages of development. 

Knowledge base is too uncertain for 
reliable modelling. 



A simple four box generic food chain model as shown in Fig. 8.1 is 
a useful screening tool that accounts for uptake from the water 
column and the lower trophic levels. This model is an improvement 
over a single BCF computation without the vast data requirement 
for a site specific food chain model. This model, however, can be 
used as a screening tool to provide a relatively rapid assessment 
for problem identification, setting priorities and planning 
further studies. A site- and species-specific model is needed in 
the event that models are to be used as regulatory tools. 

A pathway model as shown in Fig. 8.2 which includes a pelagic food 
web (phytoplankton, zooplankton, larval fish, forage fish and 
piscivorous) and a benthic food web (benthic invertebrates and 
benthic fish) may be used as part of an advanced approach to 
estimate contaminant levels in the ecosystem. 
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- The tentative diet for the generic groups shown in Table 8.1 and 
Fig. 8.2 provides insight into the type of data needed to ensure a 
proper computer model that includes the general pathways and 
coupling as illustrated in Fig. 8.2. Vhen site- and 
species-specific data become available, this generic model can be 
upgraded and calibrated for the site under consideration. 

RESEARCH NEEDS 

During the work group seminars and the formal presentations, several 
information gaps were identified as being important to the development of 
a biologically reliable model for estimating the contaminant body burdens 
in aquatic biota. 

In order to represent the impact of a proposed contaminant load on the 
biota at a particular site, it is essential that a dynamic site and 
species model be developed. At present this type of model is at the 
research stage. As a research tool the generic pathways can serve as a 
useful initial framework for the development of a calibrated and verified 
species-specific model. In order to achieve this goal there are 
information gaps that need to be filled. 

Bioavailability has been shown to be important in the uptake of 
metals. Dr. Thomann indicated that bioavailability may also be important 
in membrane transfer. Normally, dissolved contaminants are available for 
transfer across a biological membrane; however, if the contaminant is 
associated with another dissolved chemical, especially a dissolved organic 
carbon (DOC) , it may not move through the membrane of the organism as 
readily. Dr. Bierman discussed the importance of bioavailability in the 
interstitial water of sediments. More research is required to determine 
the site-specific factors that influence the bioavailability of various 
pollutants. Predictive equations are required for modelling purposes. 
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The workshop identified gut content data as one of the major 
information gaps for the Great Lakes System. This information is required 
for all potential impact areas within the Great Lakes System. From this 
information, both diet and contaminant pathways can be established. One 
of the outstanding modelling concerns that this information would help to 
answer is the relationship between the pelagic and benthic food chains. A 
considerable amount of attention in the work groups was focused on the 
applicability of size classification. More research is required to 
determine the usefulness of this concept in pelagic and/or benthic food 
chain models. Another information gap is the lack of a complete body 
burden data set for a few sites covering spatial, temporal and age 
variations of important species. This information is required for model 
calibration. 

The role of biota in the mass balance of contaminants in aquatic 
systems is not well understood. Special attention needs to be directed to 
defining the roles of macrophytes in the uptake, storage, transfer and 
release of contaminants. In addition, the importance of export due to 
emergence and harvesting on the mass balance needs to be quantified. 

Our knowledge base for the fate and uptake of contaminants with high 
K 's (log K > 7) is very limited. Dr. Mackay and his group have 
determined that these contaminants behave in a distinctive way and cannot 
be modelled adequately by the existing food chain models. This is an 
important area needing fundamental research. 

The total impact of a given discharge depends on the combined effects 
of all the contaminants in the effluent. Research into whole effluent 
modelling should be pursued. A starting point in this effort is the work 
of DiToro (1984) in which he considered modelling LC50 as a load. 
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The processes involved in the contaminant transfer from the effluent 
to the various trophic levels of the food chain are stochastic. There is 
a need to consider this stochastic nature in the modelling process. 

The accumulation of contaminants in fish depends on factors which are 
at present inadequately quantified. More research is required to 
determine the functional relationships for: growth rates of both benthic 
and piscivorous fish; respiration rates for fish; assimilation efficiency 
for food and contaminants; excretion and depuration rates for aquatic 
biota. 
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Table 8.1 FOODCHAIN DIET DISTRIBUTION IN PERCENT 

*********** * CONSUMERS ***************** # 
ORG. MACRO- PLANKTON BENTHIC ***********FISH*********** 
FOOD DET. PYTES PHYTO ZOO INV. LARVAL FORAGE PISC BOTTOM 
ORGANIC 40 4 

DETRITUS 

MACRO- 15 5 10 

PHYTES 

PHYTI- 100 15 5 10 

PLANKTON 

ZOOPLANK- 20 85 65 5 

TON 

BENTHIC 50 

INV. 

LARVAL 10 25 20 

FISH 

FORAGE 60 

FISH 

PISC. io 

FISH 

BOTTOM 
FISH 
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FIGURE CAPTIONS 

Fig. 8.1. Simplified generic food chain. 
Fig. 8.2. Suggested aquatic food chain. 
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APPENDIX A 

1. WORKSHOP AGENDA 
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3. GENERAL AND SPECIFIC QUESTIONS 



AGENDA 



AQUATIC FOOD CHAIN MODELLING WORKSHOP 



THURSDAY, JULY 9 

10:00 AM Welcome, and opening remarks (Y. Hamdy - Ontario MOE) 

10:15 AM MISA Water Quality Track (C.F. Schenk - Ontario MOE) 

10:30 AM Food Chain Modelling Framework - Purposes and 

Applications (J. A. McCorquodale - U. of Windsor) 

11:00 AM Food Chain Modelling Concepts and Requirements 
(R.V. Thomann - Manhattan College) 

12:00 NOON Lunch Break 

1:00 PM Sediment Interactions with the Biotlc and Abiotic Solids in 
the Water Column (V. Bierman - Notre Oame University) 

1:30 PM Feeding Relationships (D. Haffner/M. Weiss - U of Windsor) 

2:00 PM Uptake Kinetics (D. Mackay. - Univerlsty of Toronto) 

2:30 PM Coffee Break 

2:45 PM Work Group Orientation 

3:00 PM Work Group Sessions: i) Fisheries 

ii) Plankton 

iii) Benthos/Sediment 

iv) Macrophytes 

5:00 PM Dinner 

7:00 PM Work Group Session (cont'd) 

9:00 PM Adjourn for the day 

FRIDAY, JULY 10 

8:00 AM Meeting of Work Group Coordinators 

9:00 AM Summary and Discussion of Work Group Sessions: 

1) Fisheries 
ii) Plankton 

10:00 AM Coffee Break 

10:30 AM Summary and Discussion of Work Group Sessions (cont'd): 



Til) 
Iv) 



Benthos/Sediment 
Macrophytes/Cl adophora 



11:30 AM Overall Sumnary and Conclusions of Work Shop 

(J. A. McCorquodale - U. of Windsor) 

12:00 NOON Lunch (optional) 
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D0787B/GL-7 



July 9-10, 1987 



WORK GROU P SESSIONS 



Benthos/Sediment (Room 111) 

Coordinators: V. Bierman 
W. Scheider 



Participants : 



D. 


Andrews 


M. 


Diamond 


D. 


Mackay 


D. 


Poulton 


L, 


Wong 


B. 


Parker 


M. 


Weiss 


M. 


Kirby 



D. 


Boyd 


G. 


Johnson 


D. 


Persaud 


W. 


Wager 


M. 


Jackson 


A. 


Rodr igues 


S. 


Thornley 



Fisheries (Room 302) 

Coordinators: J. Leach 

R . Thomann 



Participants: F. Gobas 

J. MacLean 

A. Roy 

B. Shuter 
K. Clarke 



S. Klose 

D. Rokosh 

B. Keller 

G. Westlake 



Plankton (Room 314) 

Coordinators: D. Haffner 

P. Nettleton 



Participants: S. Abernathy 

P. Kauss 

N. Yan 

Z. Novak 



M. Griffiths 
A. McLarty 
E. M. Yuen 



Note : 

Macrophytes will be discussed as part of the 
benthos/sediment work group session. 

Y. Hamdy and J. A. McCorquodale will participate in each 
session to ensure continuity. 
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LIST OF MODELS USED FOR SCREENING 



A. HYDRODYNAMIC MODELS (water-column screening only) 



1. Various dilution/plume models 

2- KE MODEL - STEADY, 2"D DEPTH"AVERAGED FOR LARGE RIVERS 

3- FE MODEL " STEADY, 3'D FOR EMBAYMENTS/ LAKES 

4- RAND MODEL - UNSTEADY, 2"D FOR LAKES 

5- Spill model - lagrangian, 2-D- 



B. FATE ft TRANSPORT (hater- COLUMN, SEDIMENT, 

and target biota screening) 



1. TOXIWASP 



Primarily organics 



2 • WASTOX 



Primarily organics, with food-chain 



3. MEXAMS/MINTEQ 



Metals 
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INPUT DATA NECESSARY FOR SCREENIN6 A CONTAMINANT SOURCE 

(FOR WATER COLUMN/SEDIMENT/TARGET BIOTA CRITERIA) 



1. HYDRQDYNAMIC 

Stochastic summary for the discharge/ stage of 
the receiver* 

Current velocity measurements for calibration of 
hydrodynamic model. 



2. EFFLUENT LOADING 

BAT (EA) criteria including both concentrations and 

FLOW RATE- 



3. PARTITIONING RELATED COEFFICIENTS 

- KOH/KOC 

Organic carbon contents of sediment and target 
biota- 

4- APPROXIMATE SEDIMENT DYNAMICS 

Sediment concentrations in water column and bed 
Sediment loads from outfalls/other 

5* MEASURED AMBIENT CONCENTRATIONS 

In water column/bed/target biota, along with loading 
history/ for model calibration* 
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FATE AND TRANSPORT MODEL: TOXIWASP 
Field-acquired input data for modelling are: 

1. Particle-size distribution of suspended solids, bed sediment, and 
bed load. 

2. In-pi ace density of bed sediment. 

3. Turbidity 

4. Bed load transport rate. 

5. Pore water percolation rate. 

6. Concentrations of total contaminant and suspended solids 
throughout system (as initial conditions). 

7. Concentrations of total contaminant and suspended solids entering 
river via upstream boundary. 

8. All point-source and non-point source loading rates of total 
contaminant and suspended solids, of signficance. 

9. Hydrogen ion-hydroxide in activities (pH-pOH). 

10. Concentrations of environmental oxidants (e.g.peroxy radicals) in 
water and sediment. 

11. Organic carbon content of sediments as fraction of dry weight. 

12. Organic carbon content of the compartment biomass as a fraction of 
dry weight. 

13. Free electron activity (Eh). 

2- 2_ 

14. Concentration of major anions in water, (C0 3 , S0 4 , etc. 

2 + 2 + 

15. Concentration of major cations in water, (Ca , Mg , etc.). 

16. Alkalinity. 

17. Mineralogical composition of sediments, (Si 2 , A/ 2 2 , 
Fe 2 , 3 , MgO, CaO, Na0 2 , K 2 0, Ti0 2 , MnO, P 2 5 , etc.) 

2 

18. Concentration of Trace Constitutents, (H S, Iron, Manganese, 
PO; , F", Al 3+ , Ba 2+ , Sr 2+ , Cs + , Li + , Br", NH^ , B, 

Rb + , f , N0 3 T etc.). 

19. Concentration of Inorganic Carbon. 

20. Specific Conductivity. 

21. Bacterial population density, (both, water and sediment). 

22. Total, actively sorbing biomass. 

23. Total contaminant concentrations in specific water column and 
benthic biota, (phytoplankton, zooplankton, young-of-the-year 
fish, sport fish). 

24. Temperature (water). 

25. Wind velocity, 10 cm above water surface. 

26. Average Cloudiness, (in tenths of full sky cover). 

27. Hydrodynamic characteristics. 
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1 EXPOSURE CONCENTRATIONS 

-DISSOLVED 
-SORBED 

2 FOOD WEB . 

PREDATOR - PREY 
(BY AGE & SEASON) 

3 CONSUMPTION, 

4 ASSIMILATION 

-FOOD 
-TOXICS 

5 GROWTH RATE 

6 WET/DRY WEIGHT I LIPID 
CONTENT 

7 EXCRETION I 

8 RESPIRATION 

9 BCF 

10 BAF 

11 MOBILITY/MIGRATION/SPAWNING 

12 HARVESTING/REMOVAL 

13 DEATH & RECYCLING 

14 BIOTURBATION 

15 FEEDBACK TO EXPOSURE 
MODEL 

16 METABOLISM 
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. BENTHOS 
1. EFFECT OF BENTHOS ON PHYSICAL - CHEMICAL PROPERTIES ? 



2. EFFECT OF SEDIMENT EXPOSURES ON BENTHOS ? 



3. COUPLING OF BENTHIC AND PELAGIC FOOD CHAINS ? 



4. KEY SPECIES IN RELATIONS TO FISH - GUT CONTENTS ? 



5. ROLE OF MACROPHYTES IN FOOD CHAIN ? 
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FISH 



1. BIOLOGY : 

e.g. - PISH PRED-PREY 

- GROWTH 

- RESPIRATION 

- AS SIMULATION 

- HABITATION 

- SPAWNING 



2. CHEMICAL FACTORS 

e.g. BCF 
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PLANKTON 

1. EFFECT OF PLANKTON ON PHASE PARTITIONING ? 

- SEASONALITY ? 

- BIOMASS ? 

2. ASSUMPTIONS OF LUMPING/SPLITTING ? 
( NANO vs NET PHYTO ) 

( ZOO HERB vs CARN ) 



3. IMPORTANCE OF PHYTO BCF ON FOOD CHAIN ? 



4. IMPORTANCE OF ZOO BCF/BAF ON FOOD CHAIN ? 
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WORK GROUP 
GENERAL : 

1. EXAMINE PATHWAY SYSTEM DIAGRAM AS A GENERIC MODEL IN CONTEXT 
OP GREAT LAKES AND RIVERS. 

- COMMENT ON ITS SUITABILITY. 

2. SUGGEST AGE/SIZE CLASSIFICATION AND ASSOCIATED PATHWAYS. 



3. IDENTIFY IMPORTANT PATHWAYS AND SPECIES IN PATHWAY - WHAT 
ASSUMPTIONS ARE IMPLIED IN THESE SIMPLIFICATIONS ? 



4. COMPLETE CONSUMER FOOD MATRIX - INCLUDING REPRESENTATIVE 
DISTRIBUTION OF DIET. 



S, CHANGES NECESSARY TO UPGRADE GENERIC MODEL TO DERIVE A SPECIES 
SPECIFIC MODEL. 



6. MAJOR GAPS IN INFORMATION HINDERING FOOD CHAIN IMPLEMENTATION. 
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APPENDIX B 



WORK GROUP REPORTS 



AQUATIC FOOD CHAIN MODELLING WORKSHOP 
SUMMARY OF DISCUSSIONS OF FISHERIES WORK GROUP 



COORDINATORS: 



J. Leach 
R.V. Thomann 



PARTICIPANTS: 



F. Gobas 
J. MacLean 
A. Roy 

G. Vest lake 



S Klose 
D. Rokosh 
B. Shuter 
K. Clarke 



OPENING DISCUSSION 

The initial discussion of this workgroup centered around the issue of 
global versus specific type models of the aquatic ecosystem and the 
transfer of chemicals within that ecosystem. The participants reviewed 
the advantages and disadvantages of each direction. For global models, 
the degree of simplification is appealing since global models are 
therefore useful for screening type of activities where specific chemicals 
are to be evaluated on a broad non-site specific basis. On the other 
hand, global models do not provide specific details on individual sites 
and therefore are not useful for individual waste load allocation 
purposes. 

The present structure of the MISA program was reviewed and the issue 
of having guidelines for chemical concentrations in top predators was 
discussed. It was pointed out that some of the most significant impacts 
of toxic chemical discharges has been the closing of fisheries because of 
elevated concentrations in the top predator fish and subsequent presumed 
public health impacts if such fish are consumed. 
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APPROACH FROM TYPE OF CHEMICAL 

Following these opening discussions, the work group addressed the 
question of whether food chain modelling was required at all. As a result 
of this basic question, the group looked at dividing the chemicals under 
consideration using the octanol water partition coefficient (K ) as an 
ordering parameter. Thus the chemical population was divided into three 
broad groups: 

l0 S K o„ < 5 

log K ow = 5 to 7 

log K ow > 7 

This division of chemicals (where the bounds of such divisions are 
somewhat approximate) reflect recent work on the significance of food 
chain transfer of chemicals and subsequent accumulation up to the top 
predators. At log K < 5, food chaim accumulation is generally small. 
At log K = 5 to 7, food chain effects begin to become important and for 
log K > 7, there is uncertainty as to the actual significance of food 
chain transfers. (Summary papers presented at the workshop explored these 
concepts in greater detail.) 
Chemicals with log K < 5 

For this group of chemicals, the discussion reviewed the possibility 
of dealing with such chemicals using a simple lipid partitioning concept. 
Thus, the concentration of the chemical in a fish or any other organism on 
a lipid basis would be given by: 

' = K ow P c - CI) 

where 

v = chemical concentration in the organism (pg/kg (£p)) 

p = the fraction lipid weight, and 
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c = the water concentration (/ig/L) 

In this equation, it is assumed that K Qw is an accurate measure of the 
partitioning of the chemical into the lipids of the organism. 

The procedure then for these chemicals is therefore quite simple. An 
estimate must be made of the K Qw and the fraction lipid of the organism. 
(It was suggested after some discussion that for screening purposes, a 
conservative value of perhaps 0.20 should be used. For site specific 
.situations, the lipid fraction for individual species should be used.) 

If a guideline for the chemical concentration is available then the 
above equation is solved for the allowable water concentration, i.e. 

c = JL p (2) 

ow 

where 

v - the guideline chemical concentration in the organism 

If this concentration is greater than any other concentration 
determined from, for example, acute and chronic effects considerations on 
the ecosystem, then the guideline for the chemical concentration in the 
organisms is not controlling. If however, the water concentration 
calculated from the latter equation is less than any other water 
concentration obtained from other considerations, then the organism is 
controlling. 
Chemicals with log K > 7 

This group of chemicals was dealt with rapidly. It was concluded that 
not enough information was available to fully assess the potential for 
chemicals in this range to bioaccumulate. There is some evidence and 
theoretical justification for believing that chemicals with high K Qw 
(> 7) would not continue to bioaccumulate. 
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However, the general consensus was that this area is a major 
information gap and needs to be evaluated with field and laboratory 
research and applications efforts. If a specific request is received to 
discharge such chemicals, it was felt that the burden should be placed on 
the discharger to develop the data base necessary to justify the 
discharge. 
Chemicals with log K 5 to 7 

Since this group of chemicals has the most potential for some 
bioaccumulation, discussions centered about the appropriate model 
structures that could be used. Attempts at expanding the food chain 
structure provided to the group were not successful. It was felt that 
increasing the complexity of a generic model did serve any particular 
value although there was some disagreement in the group over this point. 
Eventually, the group concluded that it could not assign generic values to 
the various elements of a food chain model and opted in the direction of 
more site specific models. 

The basic hesitation in further expanding generic models was the 
concern that such analyses would drive the permitting process. It was 
pointed out that all of the modelling frameworks, whatever their level of 
detail were intended mainly as analyses that would provide information to 
the decision making process, but not dictate that process. 

A general consensus emerged supporting the application of existing 
food chain modelling structures to specific case studies in order to gain 
expertise in the use of such models. 

The suggestion was made also that perhaps the easiest approach was to 
calculate the chemical concentration in the top predator by simply 
applying a multiplier of 10-1000 times the concentration that would occur 



170 



from exposure to the water only. Discussion indicated that there would be 
some significant difficulty in picking the multiplier. It was also 
indicated in subsequent discussion that such a multiplier is in fact 
capable of being computed from a relatively simple generic four-level food 
chain model. 

In conclusion, the group: 

1. was split on the utility of generic models; 

2. developed a consensus on applying food chain models to site 
specific case studies; 

3. suggested several approaches for calculating chemical 
concentration in the log K Qw range of 5-7 including simple 

multiplier, or a simple four level generic model; 

4. expressed through several individuals, reservations about the 
utility of such modelling efforts in the decision making process 
and proposed instead a more specific, case by case qualitative 
and semi- quantitative determination of the importance of 
chemical transfers in the food chain. 
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AQUATIC FOOD CHAIN MODELLING WORKSHOP BENTHOS/SEDIMENT SUBGROUP 



COORDINATORS: V. Bierman 
W. Scheider 



University of Notre Dame 
Ministry of the Environment 



PARTICIPANTS: 



D. Andrews 
M. Diamond 
D. Mackay 
D, Poulton 
L. Wong 
B. Parker 
M. Veiss 
D. Boyd 
G. Johnson 
D. Persaud 
V. Wager 
A. Rodrigues 
S. Thornley 
M. Kirby 



Ministry of the Environment 
University of Toronto 
University of Toronto 
Ministry of the Environment 
University of Windsor 
Ontario Hydro 
University of Windsor 
Ministry of the Environment 
Ministry of the Environment 
Ministry of the Environment 
Ministry of the Environment 
Ministry of the Environment 
Ministry of the Environment 
Ministry of the Environment 



The Benthic/Sediment subgroup met the afternoon and evening of July 9, 
L987, charged with addressing the six general questions posed to each 
group as well as five questions specific to the subgroup. After initial 
discussions which focused on the uncertainty surrounding the use of food 
chain models, (how will they be used, under what conditions are they to be 
used) the group split into two subgroups, along the lines of expertise in 
sediment chemistry vs biology. The "sediment chemistry" subgroup was 
chaired by V. Bierman and the "biology" subgroup by W. Scheider. The two 
subgroups rejoined in the evening and the following is a summary of the 
points discussed. Most of the questions the group was charged with 
addressing were touched on. 

(i) The group consensus was that there is a coupling between benthic 
and pelagic food chains. The manner of coupling is shown in a 
generic way in Fig. B-l. Figure B-l was suggested as the basis 
for developing a generic food chain model in the Great Lakes and 
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interconnecting channels. To develop a site specific model, much 

more information would be necessary, 
(ii) The group proposed a consumer food matrix as shown in Table B-l . 

The matrix was a generic one and more information would be needed 

to apply the matrix to a specific site. The matrix will vary with 

the age of organisms (especially fish), species of organism and 

season as well as by site, 
(iii) Major gaps in formation hindering the development of a model to 

address contaminant transfer in the food chain were identified as: 

- lack of information on the role of macrophytes - to what extent 
do they take up contaminants from the sediment? 

- lack of information on food sources of biota (more gut content 
data needed) 

- lack of information on rate constants pertaining to transfer 
rates between components of food chain 

- lack of information (or perhaps a synthesis of information) on 
contaminant burdens in components of the food chain 

- lack of information on the importance of microbially mediated 
transfer processes 

- validation of any site specific model is needed before 
confidence can be placed in it. 
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Eiggre b-i 

Contaminant Movement in a Foodchain 
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Note: Contaminants dissolved in water move to all "boxes". 

Contaminants in the particulate form in water move 
to the sediment. 

Macrophytes are submergent communities 

(emergents can get contaminants from the atmosphere). 

Thi3 generic model is more applicable to organics 
than metals. 
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AQUATIC FOOD CHAIN MODELLING WORKSHOP 



COORDINATORS: D. Haffner 
P. Nettleton 



PARTICIPANTS: S. Abernathy M. Griffiths 

P. Kauss A. McLarty 

E.M. Yuen N. Yan 
Z. Novak 



GENERAL 

Plankton 

1. Effect of plankton on phase partitioning was considered to be very 

important and very dependent on the composition of the flora and fauna. 
Many plankton forms have mucilage and oil storage products, and 
considering the rapid turnover time (days) would possibly represent a key 
component in the phase partitioning of contaminants in the aquatic 
environment. 

2. Questions of lumping and splitting species, size groups, 
functions, habitats would be reduced if a size classification model were 
accepted by the MOE. 

3 k 4. As noted in 1 above, one cannot overemphasize the importance 
of phytoplankton and zooplankton on bioconcentration factors and 
luoaccumulation factors. As noted in the Lake Michigan model, this can be 
the controlling step. 

The following are the reponses for this group to the general questions 
posed in Appendix A: 

Question 1. The work group preferred to avoid using food web or food 
chain models, but recommended size class as developed by Kerr and Sheldon. 

Question 2. Size classification models as recommended above avoid 
porblems with age (e.g. species moving from one compartment to another) 
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and simplify the relationship between size class intervals. Food web 
approaches were thought to be much too complicated for any real solution. 

Question 3. Pathway and species problems are minimized using the 
approach recommended. In food webs, how do you know important pathways 
with respect to both the constancy and persistence of populations. They 
are probably quite different or at least operate on different scales. 

Question 4. Group declined to provide estimates. 

Question 5. After much discussion the group recommended research be 
oriented towards developing a hybrid of a size classification model (to 
the level of zooplankton) and food web model (forage fish to predator 
fish). This should provide managers with information on the response of 
contaminants in fish to a given control program. 

Question 6. Major gaps in information hindering food chain implements 
were : 

(a) lack of knowledge at level of species 

(b) seasonality 

(c) gut content 

(d) surface films 

(e) sediment- water flux 

(f) water- atom flux 

(g) migration 

(h) feeding habits 
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